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We analyze translationally-invariant insulators with inversion symmetry that fall outside the 
current established classification of topological insulators. These insulators exhibit no edge or surface 
modes in the energy spectrum and hence they are not edge metals when the Fermi level is in the bulk 
gap. However, they do exhibit protected modes in the entanglement spectrum localized on the cut 
between two entangled regions. Their entanglement entropy cannot be made to vanish adiabatically, 
and hence the insulators can be called topological. There is a direct connection between the inversion 
eigenvalues of the Hamiltonian band structure and the mid-gap states in the entanglement spectrum. 
The classification of protected entanglement levels is given by an integer n £ Z, which is the 
difference between the negative inversion eigenvalues at inversion symmetric points in the Brillouin 
zone, taken in sets of two. When the Hamiltonian describes a Chern insulator or a non-trivial 
time-reversal invariant topological insulator, the entirety of the entanglement spectrum exhibits 
spectral flow. If the Chern number is zero for the former, or time-reversal is broken in the latter, the 
entanglement spectrum does not have spectral flow, but, depending on the inversion eigenvalues, 
can still exhibit protected mid-gap bands similar to impurity bands in normal semiconductors. 
Although spectral flow is broken (implying the absence of real edge or surface modes in the original 
Hamiltonian), the mid-gap entanglement bands cannot be adiabatically removed, and the insulator 
is 'topological.' We analyze the linear response of these insulators and provide proofs and examples 
of when the inversion eigenvalues determine a non-trivial charge polarization, a quantum Hall effect, 
an anisotropic 3D quantum Hall effect, or a magneto-electric polarization. In one-dimension, we 
establish a link between the product of the inversion eigenvalues of all occupied bands at all inversion 
symmetric points and charge polarization. In two dimensions, we prove a link between the product 
of the inversion eigenvalues and the parity of the Chern number of the occupied bands. In three 
dimensions, we find a topological constraint on the product of the inversion eigenvalues thereby 
showing that some 3D materials are protected topological metals, we show the link between the 
inversion eigenvalues and the 3D Quantum Hall Effect, and analyze the magneto-electric polarization 
(6 vacuum) in the absence of time-reversal symmetry. 



One of the most active fields of research in recent years 
has been the study of non-trivial topological states of 
matter. The paradigm example of such a state is the 
Quantum Hall Effect, with its Integer (IQHE) and Frac- 
tional (FQHE) versions. More recently, examples of topo- 
logical phases that do not require external magnetic fields 
have been proposed, the first being Haldane's Chern In- 
sulator models Although this state has not been experi- 
mentally realized, a time-reversal invariant (TRI) version 
has been proposed and discovered! 

Recent work in the theory of topological insulators^^ 
showed that an important consideration is not only which 
symmetries the state breaks, but which symmetries must 
be preserved to ensure the stability of the state. A pe- 
riodic table classifying the topological insulators and su- 
perconductors has been created. The table organizes the 
possible topological states according to their space-time 
dimension and the symmetries that must remain pro- 
tected: time-reversal, charge conjugation, and/or chiral 
symmetriesP^ The most interesting entries in this ta- 
ble, from a practical standpoint, are the 2 and 3D TRI 
topological insulators which have been already found in 
nature!^! These are insulating states classified by a Z2 
invariant that requires an unbroken time-reversal sym- 
metry to be stable. There are several different methods 



to calculate the Z2 invariant l 3 l 7 l 9 l 12 l EH anc j a non-trivial 
value for this quantity implies the existence of an odd 
number of gapless Dirac fermion boundary states as we ll 
as a non-zero magneto-electric polarizibility in 3D J 9 * 17 * 

The current classification of the topological insulators 
covers only the time-reversal, charge conjugation, or chi- 
ral symmetries and does not exhaust the number of all 
possible topological insulators. In principle, for every dis- 
crete symmetry, there must exist topological insulating 
phases with distinct physical properties, and a topologi- 
cal number that classifies these phases and distinguishes 
them from the "trivial" ones. So far, in our discussion 
we have used the term "topological" cavalierly so be- 
fore proceeding, we should ask: what makes an insu- 
lator topological? We start by first defining a trivial 
insulator: this is the insulator that, upon slowly turn- 
ing off the hopping elements and the hybridization be- 
tween orbitals on different sites, flows adiabatically into 
the atomic limit. In most of the existent literature on 
non-interacting topological insulators, it is implicitly as- 
sumed that non-trivial topology implies the presence of 
gapless edge states in the energy spectrum of a system 
with boundaries. However, it is well known from the lit- 
erature on topological phases that such systems can the- 
oretically exist without exhibiting gapless edge modesP^ 



2 



Hence, the edge modes cannot be the only diagnostic of a 
topological phase and, consequently, the energy spectrum 
alone, with or without boundaries, is insufficient to deter- 
mine the full topological character of a state of matter. In 
the bulk of an insulator, it is a known fact that the topo- 
logical structure is encoded in the eigenstates rather than 
in the energy spectrum. As such, one can expect that 
entanglement - which only depends on the eigenstates - 
can provide additional information about the topologi- 
cal nature of the system. However we know that topo- 
logical entanglement entro py (or the sub-leading part of 
the entanglement entropy) JI2H2U the preferred quantity 
used to characterize topologically ordered phases, does 
not provide a unique classification, and, moreover, van- 
ishes for any non-interacting topological insulator, be it 
time-reversal breaking Chern Insulators or TRI topolog- 
ical insulators. However, as we will see, careful studies of 
the full entanglement spectrum^! helps in characterizing 
these statesp2H23] 

The total entanglement entropy can be continuously 
deformed to zero for trivial insulators, since the atomic 
limit to which every trivial insulator can be adiabati- 
cally continued (by the above definition) is completely 
local and has flat featureless bands. We could there- 
fore suggest that a non-trivial topological state in a non- 
interacting translationally invariant insulator should be 
defined as having an entanglement entropy which cannot 
be adiabatically tuned to zero. However, even this def- 
inition cannot be entirely correct, as the entanglement 
entropy strongly depends on the nature of the cut made 
in the system. For a single particle entanglement spec- 
trum with eigenvalues {£ } the entanglement entropy is 
determined via 

S en t =-J2(£" l0 §^ + (! - l0 §( 1 - &)) ' W 

a 

For IQH states on the sphere, the many-body wavefunc- 
tion is a single Slater determinant of occupied Landau 
orbitals, and hence an orbital cuIf^I would result in zero 
entanglement entropy since all orbitals are fully occu- 
pied or unoccupied leading to a set of {£ a } which are all 
0's or l's and do not contribute to S en t- Similarly, for a 
translationally-invariant Chern or TRI topological insu- 
lator on a lattice, a momentum space-cut would always 
give zero entanglement entropy since the Hamiltonian is 
diagonal in this basis. A spatial cut, however, would 
show mid-gap bands in the entanglement spectrum both 
the IQHE and in the topological insulator case(ie. a set 
of eigenvalues spanning the 'gap' between and 1) simi- 
lar to the ones in the real energy spectrum! 23 ! 24 ! 29 ! 30 ! 32 ! 34 ! 
These mid-gap states give large contributions to the en- 
tanglement entropy. In fact, for such states, the entan- 
glement entropy for the spatial cut cannot be made to 
vanish by any adiabatic changes in the Hamiltonian. We 
hence propose that a translationally invariant insulator 
can be classified as topological if it cannot be adiabati- 
cally connected to a state with zero entanglement entropy 
for at least one kind of cut of the system. Explicitly, an 



insulator should be characterized as topological if it has 
protected mid-gap states in the single-particle entangle- 
ment spectrum that cannot be pushed to eigenvalues 
or 1 by any adiabatic changes of the Hamiltonian. 

In the current paper we analyze the physics of insu- 
lators with inversion symmetry based on the above def- 
inition. While some of the inversion-symmetric insula- 
tors exhibit protected edge modes in the energy spectrum 
with boundaries {e.g. a Chern insulator with inversion 
symmetry), most do not. However, they can still be topo- 
logical because their entanglement spectrum for a spa- 
tial cut exhibits protected mid-gap bands of states. This 
was first pointed out for 3D strong topological insulators 
with inversion symmetry and soft time-reversal breaking 
in Ref. [55] . Although it was indicated in Ref. [55J that 
the entanglement spectrum cannot distinguish between 
a TRI and inversion invariant topological insulator, and 
one with TRI slightly broken (compared to the bulk gap) , 
we show that it can distinguish these states. In Section 
|IT] we explicitly show that inversion symmetric topologi- 
cal insulators have two types of entanglement spectrum, 
both with protected mid-gap states. The characteristic 
which distinguishes the two types of entanglement spec- 
tra (and the two cases from Ref. [2"5]) is the presence or 
absence of spectral flow. For non-trivial TRI or Chern 
insulators the entanglement spectrum exhibits spectral 
flow, very much like their energy spectrum. Heuristi- 
cally this means that the filled and empty bulk states are 
connected via an interpolating set of states localized on 
the partition between the two entangled regions. Spec- 
tral flow in the energy spectrum implies spectral flow in 
the entanglement spectrum. If time-reversal is broken for 
TRI topological insulators, or for T-breaking insulators 
with vanishing Chern number, we show that such spectral 
flow is interrupted in both the energy and entanglement 
spectra, and the occupied bands are disconnected from 
the unoccupied ones. One could then assume that, in sys- 
tems without a continuous spectral connection between 
the bulk entanglement bands, one could push all the mid- 
gap entanglement bands to entanglement eigenvalues 0, 1 
and hence to a trivial insulator with vanishing entropy 
on every cut. We find this not to be the case for inver- 
sion symmetric insulators whose set of inversion eigen- 
values changes between two inversion symmetric points 
in the Brillouin zone. In this case, while most entan- 
glement eigenvalues can be made or 1, there is a set 
of protected mid-gap states that form bands which give 
the insulator nonzero entanglement entropy (even when 
the spectral flow has been destroyed). For the case of 
time-reversal and inversion invariant topological insula- 
tors, this was shown in Ref. UHlby breaking time-reversal 
softly. We find a formula relating the number of pro- 
tected mid-gap bands in the entanglement spectrum to 
the inversion eigenvalues of the system at inversion sym- 
metric points. At inversion symmetric fc, and for cuts 
which separate the system into two equal halves, the en- 
tanglement spectrum has protected entanglement edge 
or surface modes at exactly £ = 1/2. This means that 
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there is no finite-size level repulsion (splitting) between 
these modes which is a common feature for real bound- 
ary or interface states in the energy spectrum. Even if 
the original system has other topological invariants, such 
as the Chern number, or the TRI Zi invariant which are 
all trivial, the number of entanglement edge modes can 
be nonzero. 

In Section |III| we analyze the physical response of a 
subset of these inversion symmetric insulators, and show 
important implications for the charge polarization, the 
parity of the Chern number, the 3D quantum Hall ef- 
fect, and the topological magneto-electric polarizability. 
Namely, given the set of inversion eigenvalues for the oc- 
cupied bands at all inversion symmetric points in the 
Brillouin zone, we provide explicit, compact formulas and 
complementary derivations for determining the physical 
responses which only depend on the inversion eigenval- 
ues. Additionally, as an alternative perspective, we also 
show that some of the inversion topological invariants 
are equivalent to the wavefunction monodromy, which, 
in principle is an experimentally measurable quantity. 
In particular, we show the following: in one dimension, 
the product of inversion eigenvalues over all occupied 
bands and over all inversion symmetric points is related 
to the quantized charge polarization. In two dimensions, 
the product over the inversion eigenvalues determines 
the parity of the Chern number of the occupied bands. 
In three dimensions, we show several things: first, we 
prove a topological restriction for the product of inver- 
sion eigenvalues of any insulator: it must always equal 
+1. As such, some systems are protected metals, which 
cannot be made insulating with weak scattering. Sec- 
ond, we show that, depending on the product of inver- 
sion eigenvalues in different inversion symmetric planes 
we will have 3D quantum Hall effects on different planes 
in the sample. We then show that, depending on the in- 
version eigenvalues, several inversion symmetric systems 
can exhibit a quantized magneto-electric polarizability, 
even though the Hamiltonian my not be adiabatically 
continuable to a time-reversal invariant topological insu- 
lator. 

Finally in Section |IV| we end the paper with several 
examples of such insulators and with numerical results. 
While most of the examples we chose have a topological 
response connected with an inversion topological invari- 
ant, we stress that the presence of mid-gap states in the 
entanglement spectrum is not intrinsically related to the 
presence of a non-trivial topological response. For exam- 
ple, two identical copies of a strong-topological insulator 
with inversion and time-reversal symmetry has a triv- 
ial Zi invariant and thus a trivial response. However, 
this system will exhibit protected mid-gap entanglement 
states. So, while some inversion invariant insulators have 
protected topological responses, some do not. The situa- 
tion is even more complicated: in several cases, we prove 
that inversion eigenvalues by themselves cannot uniquely 
determine the response. We can show that a nontrivial 
quantum spin Hall state and two copies of a Chern in- 



sulator each with Chern number unity have identical in- 
version eigenvalues but obviously represent very different 
states of matter. The question of relating all the inver- 
sion eigenvalues to a response function remains, in the 
cases where possible, still unsolved. 

I. PRELIMINARIES 

Let us start our discussion with some observations for 
simple two- and four-band model Hamiltonians, which 
will be referenced throughout the paper. These proper- 
ties of the models that we discuss are only dependent on 
a generic inversion symmetry and not other symmetries 
present "accidentally" in these simple models. 

2 band model. The two band model is 

H2 = £ [4 a+ V* l c *+i + h - c - + 4 (1 + m) a 3 c x ] , 

X 

where a, to are two parameters and a a are Pauli's matri- 
ces. The model is symmetric under the inversion oper- 
ation c x — >a 3 C- x . The Bloch representation of H2 takes 
the simple form: 

HzQt) = a cos k + sin(fc)<Ti + (1 + m — cos k)a 3 , (2) 

and the inversion is implemented by the operator V = a 3 : 

a 3 H(k)a 3 = H(-k). (3) 

H2 is gapped, except when m=— 2 or 0. The energy 
spectrum of the model with open boundary conditions 
is presented in Fig. [lja-b), for a=0 and to==f1, respec- 
tively. Throughout this paper, P^ will denote the pro- 
jector onto the occupied bands at momentum k, which 
is a KxK matrix (K= total number of bands), whose 
entries depend on k. 

The two special points k inv =0, tt where H2(k) is 
mapped onto itself by inversion will play a special role 
in the following discussion. We are going to examine the 
(non-zero) eigenvalues C(0) and CM of Pk=oPPk=o an d 
Pk=TrP Pk=TT , respectively, for the 3 different insulating 
phases of Hi- What we find is the following: 

1. C(0) = CM = -1, for rn > 0. 

2. C(0) = CM = +1, for m < -2. 

3. C(0) = -CM = +1, -2 < to < 0. 
We can form a Z2 topological invariant: 

xv = n w 

fci nv ,i€ occ. 

which is topologically stable since one cannot change its 
value without closing the gap of the Hamiltonian. The 
expression of \v is similar to invariants formed for time- 
reversal and inversion invariant topological insulators in 
2D and 3D.E 

For the model of Eq. [2j \v takes the values Xv=+^ 
for the insulating phases with to^[— 2,0], and xv = ~ 1 
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FIG. 1. Energy Spectra for the simple ID 2-band model 
with open boundary conditions for (a)a = m — — 1 (non- 
trivial) (b)a — m — 1 (trivial) Entanglement spectra for the 
two cases are shown in (c) and (d) respectively for a half-filled 
Fermi sea ground state with periodic boundary conditions. 



for m€[— 2,0]. A direct calculation indicates that the 
phase with x-p=— 1 displays a single end mode at each 
end of a system with open boundaries, while the phases 
with xv = +1 do not display any end modes. This sim- 
ple model indicates that, indeed, systems with inversion 
symmetry do posses non-trivial topological phases that, 
apparently, can be classified by the xv invariant. As we 
shall see in the following sections, xv can be- linked to a 
physical response of the system but it cannot completely 
classify the topological phases of a system with inversion 
symmetry, even in ID. 

4 band model. It is instructive to repeat a similar anal- 
ysis on a 4-band inversion symmetric model. For this we 
use the following Hamiltonian, written directly in the 
Bloch representation: 



H 4 (k) = sin(fc)f i + sin(fc)f 2 
-(1 - m - cosfc)f + 5T 2 4 +ecos(fe)(l + f ) 



(5) 



where Ti = <r z <g> t x , T 2 = 1 <S> r y , T = 1 ® t z , and 
F24 = o~ x (E) t z . The Pauli's matrices T a ,a a act in the 
orbital and spin spaces, respectively. H 4 (k) is symmet- 
ric under inversion, which is implemented by V = Tq, 
and is gapped except for a few values of the parameters 
8, e, and m. Note that this system also has an acciden- 
tal time-reversal symmetry with T = (io~ v <g> 1)K, but 
this can be broken without affecting the stability of the 
topological state or removing the mid-gap modes in the 
entanglement spectrum. The two lower energy bands are 
assumed occupied, and in this case PqVPo and P^VP^ 
are 4x4 matrices, each displaying two non-zero eigenval- 
ues Ci(0) and Q(tt), i—1,2. 



We are going to present the inversion eigenvalues for 
the insulating phases of the model. There are six such 
phases (we discuss only five of them) and their energy 
spectra with open boundary conditions are shown in 
Fig. |2ja-e). Choosing representative values for the pa- 
rameters, we find: 

Case 1) b = 8 = e = and M < 0: 



Ci(0) = Ca(0) = +l 
CiW = CaW = +1. 



(6) 



and consequently xv = +1- 
Case 2) b = 8 = e = and M 



0.5: 



Ci(0) 



■■ C 2 (0) 
C2M 



(7) 



and consequently xv = +1- 

Case 3) b = 1, 5 = 0.7, e = and M =< 0: 

Ci(o) = -i, c 2 (o) = +i 

GM = CaW = +1. 

and consequently xv = — 1- 

Case 4) b = 1,6 = 1.7, e = and M = 0.5: 

Ci(0) = -1, Ca(0) = +1 
Ci(7r) = +1, CaW = -l. 

and consequently xv = +1. 

Case 5) b = 1, 5 = 1.7, e = 0.7 and M = 0.5: 

Ci(0) = +1, Ca(0) = -1 
Ci(7r) = +1, CaW = -l. 



(8) 



(9) 



(10) 



and consequently xv = +1- 

The 4-band model reveals a far richer internal struc- 
ture. Case 1) can be identified with a trivial topological 
phase and, based on the value of xv and on the pres- 
ence of end modes seen in Fig. [2jc), one will be inclined 
to classify case 3) as a non-trivial topological insulator. 
But one will have clear difficulties with labeling cases 2), 
4) and 5). This is a clear indicative that xv alone is 
not enough for a full classification and that additional 
topological invariants are needed for a complete picture. 

It is instructive to consider the atomic limit of the 
model. By the atomic limit we mean the limit of the 
adiabatic process in which the hopping terms between 
different sites are tuned to zero. Since the bands are 
dispersion-less and completely local (disentangled) in this 
limit, it makes sense to talk about the parity of an en- 
tire band (or orbital), since its inversion eigenvalues at 
k = 0, 7T are identical. For a model with 2 occupied 
bands, the atomic limit can lead to the following cases, 
depending of how the occupied atomic orbitals behave 
under inversion: two occupied bands of parity + (labeled 
++), two occupied bands of parity — (labeled ), and 



5 




(f) 



0.5 



(g) 

0.5 



(h) 

5 1 : 

0.5 



(i) 

0.5 




FIG. 2. (a,b,c,d,e) Energy spectra for Hi in cases 1,2,3,4,5 respectively with open boundary conditions. (f,g,h,i,j) Entanglement 
spectra for H4 in cases 1,2,3,4,5 respectively with periodic boundary conditions. 



one band of parity + and one band of parity — (labeled 
H — ). These three options give the complete classifica- 
tion of the trivial inversion symmetric insulators with 2 
occupied bands in ID. Now, a direct calculation will show 
that Cases 1), 4) and 5) can be connected to their atomic 
limits without closing the insulating gap and that Case 
1) can be identified with H — h trivial insulator, while both 
Cases 4) and 5) can be identified with the H — trivial in- 
sulator. Note that Cases 4 and 5 can be adiabatically 
connected to each other without closing the bulk insulat- 
ing gap. The trivial insulator occurs in our 4-band 

model if we take the large m limit. 

Based on the above discussion and on the ab- 
sence/presence of the end modes in Fig. [2j we can con- 
sider the Cases 1, 4, 5 as completely trivial and Case 3 as 
non-trivial, but Case 2 is still uncharacterized. It cannot 
be continued to the trivial atomic limit without closing 
the bulk gap, it displays end modes, yet \v = +1- To dis- 
tinguish this phase we must carefully consider the parity 
eigenvalues. We see that when the fc inv points are consid- 
ered separately, Xv(k nv ) = Y[ ieocc . d(k nv ) = 1- Thus, at 
each fc inv there are an even number of bands with nega- 
tive parity eigenvalues. In this situation, when for each 
fc inv the local product over the eigenvalues of the occupied 
bands is trivial (+1) we can define a second invariant 

x^ 2) ^n fcl „ v , eocc ./ 2 c I (fc 1 n V ), (11) 

where the product over bands is defined to be the product 

of half of the bands with negative parity eigenvalues at 

each fc inv . Note, that we do not require there to be an even 

number of filled bands, just an even number of negative 

parity eigenvalues. Out of the five cases discussed above, 
(2) 

the x-p invariant can only be defined for cases 1 and 2 
(2) 

for which x-p is trivial/non-trivial, respectively. As we 



shall see, the x-p invariant is relevant and important for 
inversion symmetric insulators in 3D. 

The invariant x—\ n i~ n 2\-> where n\ and n 2 are the 
number of negative parity eigenvalues at k=0 and 7r, re- 
spectively, is also useful for classifying the inversion sym- 
metric insulators. This invariant generically indicates 
how many times the insulating gap must close when one 
takes the atomic limit. As we shall see, x also gives the 
number of robust mid-gap modes in the entanglement 
spectrum localized on a single cut-boundary. For our 4- 
band model, x=0 for Cases 1), 4) and 5); %=1 for Case 
3) and x=2 for Case 2). 

To conclude, the 2- and 4-band explicit models show 
that the systems with inversion symmetry can display 
topologically distinct phases, i.e. they cannot be con- 
tinuously deformed into one another without closing the 

insulating gap. These phases can be characterized by the 

(2) 

different values of %, xp or Xp (when they can be con- 
sistently defined). In the remainder of the manuscript, 
we will try to elucidate the physical signatures of these 
topological phases and to give the physical meaning of 
the invariants mentioned above. For this, we will investi- 
gate the properties of the entanglement spectrum, which 
we demonstrate to contain robust edge modes, and the 
physical responses of the inversion symmetric insulators. 



II. ENTANGLEMENT SPECTRUM OF 
TOPOLOGICAL INSULATORS WITH 
INVERSION SYMMETRY 

In this section we discuss the bi-partite single-particle 
entanglement spectra for inversion symmetric topolog- 
ical insulators. Previous work on entanglement spec- 
tra in translationally invariant topological insulators was 
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carried out in Refs. [531 HH 133 and [MJ where it was 
shown that the primary contributions to entanglement 
arise from states localized near the spatial cut between re- 
gions A and B. Additionally the entanglement spectrum 
of disordered Chern insulators has been investigated in 
Ref. |32l The first indication that the presence of in- 
version symmetry is important for the structure of the 
entanglement spectrum was presented in Ref. 1291 Here 
it was shown that, while the physical edge spectrum of 
a time-reversal and inversion invariant topological insu- 
lator is gapped in the presence of an added Zeeman field 
(which does not close the bulk gap), the entanglement 
spectrum still contains a gapless mode. The authors of 
that work link the existence of mid-gap states for each 
cut in the entanglement spectrum with the existence of 
a 6 = 7r vacuum characteristic of a TRI non-trivial topo- 
logical insulatorP Although, as we mentioned before (and 
will discuss more in Sec. Ill), there is not always a di- 



rect and unique connection between the physical response 
and protected states in the entanglement spectrum, this 
was an essential indication that inversion symmetry could 
support topological states and that the properties of the 
entanglement spectra were closely connected with inver- 
sion symmetry. 

We start this section by detailing how to obtain the en- 
tanglement spectrum for noninteracting insulators. We 
then look at the entanglement spectrum of topological in- 
sulators and show that there arc two fundamental proper- 
ties which may be present (i) protected mid-gap states at 
entanglement eigenvalue £ = 1/2 and (ii) spectral flow in 
the entanglement spectrum. In the presence of inversion 
symmetry, there can be mid-gap states in the entangle- 
ment spectrum, and these may or may not be connected 
to the entanglement bulk band edges via a spectral flow 
pattern. Hence there are two distinct types of non-trivial 
entanglement spectra. One example that we will see is 
TRI topological insulator parent states which are inver- 
sion symmetric but which may have time-reversal slightly 
broken. The time-reversal invariant case has both pro- 
tected mid-gap modes and spectral flow, while the T- 
broken case only has protected mid-gap modes. 

For now we focus solely on insulators with a generic 
inversion symmetry and show its consequences on the 
entanglement spectrum. First, we show that if the system 
is cut exactly in half, then there can be mid-gap states in 
the entanglement spectrum located exactly at a value of 
1/2. This is equivalent to the statement that the mid-gap 
eigenvalues of the flat band Hamiltonian, for cuts exactly 
in half, exhibit no finite-size level repulsion. We then 
give an expression for the number of 1/2 eigenvalues in 
the entanglement spectrum as a function of the numbers 
of negative inversion eigenvalues at inversion symmetric 
points. Ref 29 also points out the existence of multiple 
exact midgap states in the entanglement spectrum but 
does not relate them to the difference of the inversion 
eigenvalues between inversion symmetric points. 

In everything presented below, it is very important to 
clarify that by a spatial cut in a translationally invariant 



system we mean a cut between primitive unit cells. This 
point is important when we consider systems with par- 
tially broken translation symmetry e.g. the dimcrizcd 
models in Sec. |IV| The physics is independent of the 
choice of the unit cell. For example, in a multi-orbital 
system with a one-site unit cell, the cut should not be 
made through the orbitals on the same site. 



A. Obtaining the Entanglement Spectrum of an 
Insulator 



All of the models we study are free fermion Hamil- 
tonians. To find the single-particle entanglement spec- 
trum we use Peschel's method.^ We begin by assuming 
a quadratic Hamiltonian of an insulator with a = 1...K 
quantum states per site which is translationally invariant: 



h = J2 c Lk H M k ) c P,k 



(12) 



and the canonical transformation U that diagonalizes it: 

U^HU = Diag(£ n ). (13) 

U is the matrix of eigenvectors u n (k) of energy E n (k): 

U(k) = (u 1 (k),u 2 (k)....u K (k)) (14) 

where each u n {k) is a if-component vector. In general, 
we will use k to denote the wave-vector of components 
kx-, ky, etc.. 

The relationship between the normal mode operators 
"/pk and the electron creation operators is 



Cafe = U a p(k)jp k = U™(fc)7„ fc 



(15) 



To calculate the single-particle entanglement spectrum 
we simply need the correlation function: 



C 



a/3 



(16) 



where c\ a creates an electron in state a at site i. The 
expectation value is taken in the ground state. We can 
view this correlator as a matrix Cy, with entries that 
depend on i and j. We have: 

n a P — \^ Jkii-ik 2 j /J rai \ 
fci ,fe 2 



k n£occ. 

= Y^e lk ^-^Pf 1 (17) 



or more compactly 



(18) 



We want to make a translationally invariant cut along the 
y-direction so that k y is still a good quantum number (k y 
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is a short-hand notation for all the momenta parallel to 
the cut, so we are implicitly also treating systems in 2, 3 
dimensions). Thus we have: 

where L is the total number of sites along the cut. Follow- 
ing Peschel, for the entanglement spectrum, we restrict 
i,j to be in region A, which is an explicit cut in position 
space. There are several physical choices for cuts, but for 
topological insulators we will show that a spatial cut can 
distinguish between topological and trivial insulators. 

As an aside, note that for an insulator the spec- 
trally flattened Hamiltonian matrix where the states 
above/below the gap are flattened to energies +1/2 and 
— 1/2 respectively is given by: 

H &&t {k x ,k y ) = \-P K , kv . (20) 

The two above expressions of the correlation function 
and of the flat band Hamiltonian explicitly show that 
the entanglement spectrum, i.e. the eigenvalues of the 
restricted C°"? are identical to the energy levels of the 
flat band Hamiltonian with open boundaries in region A 
shifted by a constant. As such, if the flat band Hamilto- 
nian is topological (i.e. has protected edge states), then 
immediately we know the entanglement spectrum will 
have states localized on the cut. This agrees with the 
results of Ref. [30l The interesting thing is that the flat- 
tened Hamiltonian can have mid-gap states even when 
the un- flattened one does not. 



B. Properties of the Entanglement Spectrum 

We would like to first get an intuitive idea of how the 
entanglement spectrum of an insulator should look. The 
one-body correlation function over the full system (not 
only over region A) is a projector. It is the real space 
representation of the projector onto the occupied bands, 
and as such only has the eigenvalues and 1. This is 
shown explicitly in Appendix [XJ When we make a cut, 
the eigenvalues of the one-body correlator deviate from 
0, 1 but most of them only deviate slightly. However, 
in the topologically non-trivial case, we must get entan- 
glement "edge" -modes similar to the edge states, but 
localized on the entanglement cut, because we are re- 
ally diagonalizing the spectrum of the open boundary 
flat band Hamiltonian of a topological insulator. 

From now on we choose to cut the system exactly in 
half. A cut exactly in half will enable us to show the ex- 
istence of exact degeneracies rather than exponential de- 
generacies only arising in the thermodynamic limit. Our 
choice of cut does not matter in the thermodynamic limit, 
where it cannot physically matter whether we make a cut 
exactly in the middle or away from the middle. It is how- 
ever the case that if we cut the system in two identical 



halves, we can prove things exactly, otherwise we can just 
give arguments. 

The one-body matrix, computed in the basis <4 Q |0), 
takes the block form 

c < 21 » 

where Cl is the matrix of the left-half (the one we di- 
agonalize for the entanglement spectrum) (C^ij = Cy , 
i, j G A; Cr is the matrix of the right-half (Cr)^ — C I "' 3 , 
i,j € B] C LR is the left-right coupling, (C LR )ij = C"f , 

i E A, j e B; with C RL = C\ R . Since = C i+n>j+n , 
the following extra property is true if the cut is exactly 
symmetric (in which case a proper translation of A gives 
B): 

C R = C L . (22) 

We also find that Clr = C LR . Moreover, the projector 
property C 2 = C: 

( c L C LR \(C L C LR \_(C L c lr \ 
Wlh Cl ){Crl C l )-{c{ r C l ) (23) 
gives the following additional identities: 

Cx(l — C L ) = C lr Cl R , 

ClrC\ r = C\ r C L r (24) 

ClClr + ClrPl — Clr- 

Using the last equation, if t/j is an eigenstate of the en- 
tanglement spectrum matrix Cl with eigenvalue (proba- 
bility) p: 

C L ^ = pip (25) 
then ClriP is also an eigenstate with eigenvalue 1 — p: 
C'lClr^ = Clr^ - ClrCl^P 

(26) 

= c L Rip - pC L Ri/j = (i - p)c L Ri>- 

lip = 1/2, ifi and ClriP have the same 1/2 entanglement 
probability, but as we shall see, this does not automati- 
cally mean that the p = 1/2 entanglement probability is 
doubly degenerate because ip and Clr4> are n °t linearly 
independent, in general. 

1. Properties of Entanglement spectrum with time-reversal 
symmetry 

The entanglement spectrum maintains the symme- 
tries of the original Hamiltonian. For example, for 
time-reversal symmetry of the original Hamiltonian 
TH{k)T- 1 = H(-k) (equivalently, TP k T- x = P_ fc ): 

TCi^T- 1 = Y. kx Te^-VP^T- 1 = 

= E fcx e- it «(H)rA„ t „r 1 = 

= £ fex e-M^)P_ fex ,_ fey = 
= E fc . ^- j) PK,-k y = C l3 {-k y ) (27) 



so we see that the correlator also has time-reversal sym- 
metry, and for spin 1/2 particles for which T 2 =— 1, 
the entanglement levels come in pairs at fc and —k. 
Thus, there are entanglement Kramers' doublets at time- 
reversal invariant points where fc = — fc mod G where G 
is a reciprocal lattice vector. 



C. Inversion Symmetric Topological Insulators 

In this section we give explicit arguments that the en- 
tanglement spectrum of an insulator with inversion sym- 
metry (and without any other symmetry) can have mid- 
gap states pinned at exactly 1/2, without level repulsion 
when cut exactly in half. An integer number of such 
modes is robust without splitting, so the classification 
of the entanglement spectra of insulators with inversion 
symmetry is given by an integer Z (compare with the Zi 
case where an even number of modes would be unstable) . 
As an example, in ID, if a bulk insulator has a number 
m of filled bands with negative inversion eigenvalues at 
k = and a number ni at k = w, we give explicit ar- 
guments that the entanglement spectrum for a system 
with periodic boundary conditions (when the system is 
cut exactly in half) will have 2\ni — H2I protected mid- 
gap modes at exactly 1/2. In more than ID, there will 
be conserved momenta parallel to the cut (say k y ), for 
an insulator cut in the x direction. When cut exactly in 
half, there will be 2|m — n%\ zero modes situated at the 
Ky — —Ky mod G y for which, in the periodic bulk (be- 
fore the cut), there were ni negative inversion eigenvalues 
at (k x , k y ) = (0, Ky) and 712 negative inversion eigenval- 
ues at [kxi ky ) = (n,Ky). We illustrate this explicitly 
with several examples in Sec. IV 



1. Properties of Entanglement Spectrum With Inversion 
Symmetry 

With inversion symmetry: 

VH(k)V- x = H(-k), V 2 = l; V = V~ Y (28) 

we can define a unitary matrix Hy(fc) which connects the 
bands at k and — k: 



\u i (-k))=B*Jk)V\u j (k)) 



(29) 



where the indices i,j run over the occupied bands 
1,...,N. In fact, by performing simple band crossings 
between the N bands below the gap (which does not in- 
fluence the physics in the gap which depends only on the 
ground state) , we can make the bands non-degenerate in 
which case we can use B*j(k) = e*+W5ij, but we do not 

need to choose this gauge here. Since i\ = r/(EF — H(k)), 
where rj{x) is the Heaviside function, we have: 



VP k V = P-k, 



(30) 



which can be used to show: 

P&ijP = e tk(l ~ 3) P-k = Cji = C\ r (31) 

We now want to relate the appearance of these 1/2 
eigenvalues with the inversion eigenvalues of the occu- 
pied bands. We first consider the one-dimensional case 
where we will be able to infer the behavior of the insu- 
lator just from k x = 0, n. In principle only two sites in 
the x direction should be enough to reveal the physics. 
Of course, with just two sites, our cut has to be made 
right in the middle of the two-site problem, i.e. we are 
computing the entanglement spectrum of one site vs. the 
other site. This seems a bit problematic at first because if 
we are looking for the properties of the energy spectrum 
in a topological insulator phase the wavefunctions of the 
states localized on each end will overlap and the degener- 
acy of these low-energy end states will be lifted because 
of the small size. Crucially, we show that the flat-band 
Hamiltonian does not have such finite-size eigenvalue re- 
pulsion between the edge modes even when these modes 
rest on top of each other on the same site. That is, even 
if we bring the ends close to each other e.g. on the same 
site (which is the meaning of the one-site entanglement 
spectrum), it is still true that the end modes do not ex- 
hibit level repulsion and are degenerate. This statement 
is true in higher dimensions where the end states be- 
come propagating edge and surface states. We prove this 
statement for several particular cases, which indicate it is 
true in the thermodynamic limit. We first show that for 
one occupied band (we do not particularize to a specific 
model) , there are two mid-gap modes at exactly 1 /2 if the 
inversion eigenvalue at k = 0, 7r are opposite. Then we 
repeat this procedure for a chain of four sites cut in half. 
We then show that for two occupied bands (we do not 
particularize to a specific model), there are two mid-gap 
modes at exactly 1 /2 if there is one inversion eigenvalues 
at k = 0, 7r opposite (while the other two are the same), 
whereas there are four mid-gap modes at exactly 1/2 if 
both inversion eigenvalues at fc = are opposite from the 
ones at fc = n (i.e. at one momentum both are negative 
and at the other momentum both are positive) . We again 
do this for a chain of two sites cut in half, then for a chain 
of four sites cut in half. The main conclusion to be drawn 
from this is that in the flat-band Hamiltonian (entangle- 
ment spectrum), these mid-gap modes do not experience 
eigenvalue repulsion. It is physically clear that, although 
our proofs are only for two and four site flat band Hamil- 
tonians (cut in half for the entanglement spectrum) , level 
repulsion will not set in for larger systems: level repulsion 
between edge modes gets weaker as the distance between 
them is increased. Finally, at the end, we look at the gen- 
eral case of N occupied bands for the two-site problem 
and prove that the number of 1 /2 modes in the entangle- 
ment spectrum is 2 1 n-i — ri2 1 where n\, 77,2 are the number 
of negative eigenvalues at k = 0, tt. As there is no level 
repulsion when all modes are spatially on top of each 
other, we do not expect level repulsion when the num- 
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ber of sites is increased to the thermodynamic limit. We 
check this numerically for several examples with larger 
system sizes (e.g. 100 sites). Our exercise shows that 
time-reversal invariant insulators with inversion symme- 
try (or even the case with T slightly broken) are not the 
only inversion symmetric topological insulators with pro- 
tected entanglement mid-gap states. These are but one 
of a whole series of inversion symmetric insulators with 
mid-gap entanglement modes. 



2. One Occupied Band, Two-Site Problem 

First we look at a generic case with one occupied band, 
two-sites, and periodic boundary conditions. In this case, 
/c-space contains only the points k—0, n. The wavefunc- 
tion of the occupied band is \ipi(k)}: 

H(k)ipi(k) — e(k)tpi(k) (32) 
with inversion eigenvalues: 

v \i>x (o)) = c(o) |Vi(Q)> , v |^(tt)) = CM |Vi W) . 

Since V is a unitary operator which squares to unity, we 
have: =V(= V~ v ) and by taking scalar products in 
the above we have: 



(C(o)-cW)(^i(o)|ViW) = o 



(34) 



Hence if £(0) = — ((ir) (the eigenvalues can never be zero 
due to dctV = 1) we have (V>i(0)|^i(tt)) = 0. Notice 
that the Hamiltonian H(k) does not impose any restric- 
tions on the wavefunctions at different momenta k, i.e. 
at fc=0 and k=n we are effectively diagonalizing inde- 
pendent Hamiltonians. What allows us to relate wave- 
functions at k = 0, 7r is that they are both eigenstates of 
the same matrix V (it is important to recall that V is 
k- independent) . 

For the two-site problem (i = 1,2): 



(35) 



The eigenstates of the original Hamiltonian have opposite 
inversion eigenvalues then per the above: 



PolVlW) = Ar HMo)> =0 



(36) 



which means that ^i(0) and ipiiyi), the original Hamil- 
tonian eigenstates, are also the eigenstates of the entan- 
glement spectrum, with two eigenvalues at 1/2: 

C L |Vi (0)) = \ |^i(0)) ; C L \M*)) = \ l^i W) • (37) 

We see that the original Hamiltonian can change, leading 
to a change of ip\(k), but as long as the inversion eigen- 
values remain fixed and opposite to each other, and as 
long as we can take the flat-band limit (both of which 
mean no gap closing), the eigenvalues of the entangle- 
ment spectrum will be fixed at 1/2. It does not matter 



what the actual explicit model for H(k) is. If the inver- 
sion eigenvalues at k — 0, 7r are not opposite, there is no 
reason why (■0 1 (O) 1-01 C 71 ")) = 0, and the 1/2 modes might 
not exist or will not be protected. 



3. One Occupied Band, Four-Site Problem 

With one occupied band with opposite inversion eigen- 
values, the two-site Hamiltonian has exact 1/2 modes in 
the entanglement spectrum. This is the first indication 
that the modes are stable and experience zero level re- 
pulsion. We now show that the generic four-site Hamil- 
tonian, with one occupied band also has exact 1/2 modes 
without level repulsion. This strongly suggests that these 
modes are stable in the thermodynamic limit, as long as 
the entanglement spectrum is computed for a system cut 
exactly in half. But, in the thermodynamic limit we know 
there can be no physical difference between a cut in half 
and any other cut except for exponentially suppressed 
finite-size level splittings. Thus there will be asymptotic 
zero modes in the thermodynamic limit regardless of the 
cut. 

For the four-site problem, fc-space contains four mo- 



menta kj 



JJL 
2 



J 



0,1,2,3. Call the occupied eigen- 



state of the Hamiltonian, as above \ipi(k)), H(k)ipi(k) = 
s(k)ipi(k). For a system cut in half, the entanglement 
spectrum is given by diagonalizing the matrix Cy = Cj_j 
with i, j — 1, 2: 



r , _ l C\i C12 
u i - I At n 



(38) 



We have: 

C\\ = C*22 
Cl2 = 



-■[P Q +P*+P«+Ps«) = d a (39) 



-(P -P w + iP% -iP*. 



Ci. 



(40) 



The eigenstate of Cl corresponding to the eigenvalue £ 
takes the form (i^AtiPb), which satisfy the equation: 

CoVa + Ci^b = CiVm + C Q i{>B = (41) 

Due to the presence of inversion symmetry, irrespec- 
tive of the inversion eigenvalues, we showed before that 
VCijV = Cji — Cjj, which renders the second equation 
ofEq.[4TJ 

C x V^ a + C ViI>b = £,V^b- (42) 
We see that this is consistent with the first equation of 



Eq. 



41 



if ipB = TriPipAi with m 



equation to solve is then: 

{C + mCxV)^A = 



1. The eigenvalue 
(43) 
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Also because of inversion symmetry, we have that: 

P^^VPkV. (44) 
For the one-band problem, we know that: 

P r = c(o)p ; P*v = (45) 

where C(0), CM are the inversion eigenvalues at k = 0, 7r 
of the occupied band ip\{k). Hence to find the entan- 
glement spectrum we need to diagonalize the following 
operator: 

F=i[(l + mC(0))P + (l-mC(7r))A r 

+ P^V(V + im.) +VP^(V -im)]. (46) 
For the half-mode, we pick an ansatz: 

il>A = ai/> 1 (0) + tyi(*)> ( 47 ) 

which we show can diagonalize F for an appropriate 
choice of a, b: 

a=-(l + imC(7r))(Vi(f)|ViW) 
6=(l + imC(0))<Vi(f)|Vi(0)). 

This choice of a and b makes (Co+mCiV)tpA independent 
of both Vi(f ), and -PVi(f ) in general (i.e. [PkV(V + 

im) + VP^CP — im)]ipA — 0), as it should in order for 
our ansatz to be an eigenstate. With this choice of a, b 
we find (by taking ±[(1 + mC(0))P + (1 - mC^))P^ A ) 
that in general the entanglement spectrum eigenvalue is 
dependent on (-i/'i(O)|'0i(7r)), and hence the mode is not 
fixed at 1/2. However, if the inversion eigenvalues at k — 
0,7r are opposite ((n) = — ((0) then (-0i(O) (7r)) = 
and we find the eigenvalue of the entanglement spectrum 
of our ansatz to be: 

£=I(l + mC (0)) (49) 

Recall that have the liberty to choose the values of m = 
±1, which is equivalent to saying to — ±£(0). If we 
pick to — C(0), then our ansatz gives an eigenstate with 
entanglement eigenvalue equal to exactly 1/2. The other 
choice leads to £ = 0, so is of no interest to us. The 
eigenstate at 1/2 is (^a,^oP^Pa), where: 

yjA = i(Ml)\M*))M0) + (Ml)\M0))M*) (50) 

In the pathological case when (V>i(f )IV'i( 7r )) = 
(V>i(f )|Vi (0)) = 0, both Vi(0), Vi(tt) are 1/2 modes, but 
in general only their combination ipA gives a robust 1/2 
mode. 

A second stable 1/2 eigenvalue can be found by picking 
to = ((it) = — C(0), but using a different ansatz: 

^=aVi(f)+Wi(f). (51) 



With this choice for to, the F matrix to be diagonalized 
is: 

\[P^V{V iC(0)) + VP^V + iC(0))]. (52) 

After straightforward calculations, we find the exact 1/2 
mode to be: 

^ = Vi(f)+<C(0Wi(D- (53) 

Thus, for a completely generic Hamiltonian and its eigen- 
states we have shown that for one occupied band, if the 
inversion eigenvalues at 0, 7r in the bulk are the oppo- 
site of each other, there are two exact midgap states at 
1/2 in the entanglement spectrum. These mid-gap states 
have linearly independent eigenstates: (tpA, G(0)VipA), 
(V'ai ~ C(O)PV'a) an d V'A, iP'a as above. The standard in- 
tuition about interface or boundary states is that the 
eigenvalues repel less as the length of the system is in- 
creased, and the entanglement spectrum for the case an- 
alyzed here will have two exact 1/2 modes in the ther- 
modynamic limit, when cut exactly in half. Numerical 
simulations on specific models agree. In the thermody- 
namic limit, there can be no physical difference between 
a half cut and a cut away from half, and the levels will 
asymptote to 1 /2 in the infinite size limit if the cut is not 
exactly in half. 

4- Two Occupied bands, Two-Site Problem 

On our way to the most general case we now analyze 
the two-site, two occupied band problem. This follows in 
the same fashion as the previous example, except that we 
now have an extra complication. Namely, there are more 
options for the sets of occupied-band inversion eigenval- 
ues. 

We have the two occupied bands of the original Hamil- 
tonian: 

H(k)Mk) = e!(k)Mk);H{k)Mk) - e 2 (k)Mk) (54) 

For k = 0, 7r we generically have: 

(Vi(0)|^(0)> = <ViWI^W)=0. (55) 

We denote the inversion eigenvalues for -01(0), "02(0), 
"01 (tt), 02 (tt) by Ci(0),C2(0),Ci(7r),C2(7r) respectively. 
The order of the occupied inversion eigenvalues at each 
inversion symmetric momentum can be changed without 
affecting the topological structure so we assume that all 
negative inversion eigenvalues are listed first. We now 
calculate the number of midgap 1/2 eigenvalues in the 
entanglement spectrum. If the number of negative in- 
version eigenvalues at k = is the same as the number 
of negative inversion eigenvalues at ir (which means the 
number of positive inversion eigenvalues is also the same) , 
it is easy to prove that in general there are no protected 
1 /2 modes because the entanglement eigenvalues depend 



11 



on the overlap of bands at the two inversion symmetric 
momenta (cf. Eq 61 1. If the number of negative eigenval- 



ues at the inversion symmetric points is different, then 
we distinguish two cases: 

Case 1: The number of negative eigenvalues at k = 
differs from the number of negative eigenvalues at k = tt 
by ±2 (i.e. they are both different): 

Ci(o)CiW = Ci(o)C 2 W = C 2 (o)CiW = C 2 (o)C 2 M = -1 

(56) 

implies: 

(^l(0)|Vl(7T)) - (V 1 (0)|^ 2 (7r))=0 

(ih(P)\M*)) = (Mo)\Mn)) = o. (57) 

The one-site entanglement spectrum obtained by cut- 
ting the system in half is obtained by diagonaliz- 
ing the operator Ct, = (Po + P w )/2 where P = 

ELi \Mo)) M0)\, h = Eli (V'iWI- Due to 

their inversion eigenvalues, eigenstates at tt have zero 
eigenvalue under the projector at (and vice-versa) but 
unit eigenvalue under the projector at tt. We see that the 
modes at 1 /2 in the entanglement spectrum are given by 
exactly the occupied eigenstates of the original Hamilto- 
nian ipi(0), "0 2 (O), V'iWj ^ 2 (tt). There are exactly 4 of 
them, twice the difference between negative and positive 
eigenvalues at the two inversion symmetric points. 

Case 2 The number of negative eigenvalues at k = 
differs from the number of negative eigenvalues at k = tt 
by ±1: this implies that at one k point, both inversion 
eigenvalues are identical. Without loss of generality, let 
this point be k = tt and let the eigenvalue products be: 

Ci(0)CiW = Ci(0KaW = -l 

C 2 (o)CiW = C 2 (o)C 2 W = i (58) 

which renders the following inner products to be zero: 

(M0)\M0)) = Wi(o)|^W) = (M0)\Mk)) = o. 



Consider the eigenvalue problem: 



(59) 



We expand the state ipA into the (non-orthogonal) set of 
eigenstates ipi(0), ^ 2 (G), ipi(ir), ifa(ir): 

\M = oi |Vi(0)> + a 2 |^(0)) + h \if!i(ir)) + b 2 \M*)) 

(60) 

As ^i(O) is orthogonal with all the other eigenstates at 
both k = 0, tt since it has a different inversion eigen- 
value, it is then obvious to see that the first 1/2 mode 
solution is (01,02,61,62) = (1,0,0,0). To find another 
1/2 mode, we must expand in the three remaining eigen- 
states: \iP A ) = a 2 |V> 2 (0)) + 61 |-0x (tt))+ 62^2 (tt))- There 
is a slight complication with this expansion since noth- 
ing guarantees that the states 1^2 (0)) , IViC 71 ")) > IV' 2 ( 7r )) 
are orthogonal: in fact, in the generic case, they are not. 
Moreover, it is not clear that they are even linearly in- 
dependent. We will assume that the states are linearly 



independent. This is a perfectly valid procedure since if 
the |^2(0)) , \ipi(n)) , ^(ct")) are not independent, we will 
simply get a non-trivial nullspace. However, the non-zero 
eigenvalues are still good eigenvalues of the entanglement 
matrix. The matrix to diagonalize is: 

\ (v 2 (o)|Vi«) iMVlM*))' 







(61) 

with an obvious 1/2 eigenvalue for the state (a 2 , 61, 62) = 
(0,-<^2(0)|V 2 W),<^ 2 (0)|Vi(7r))) (01 =0). In the non- 
generic case when (ip 2 WlV'aW) = (^ , 2 (0)|V'i( 7r )) = 0, 
■02(0) is the other 1/2 eigenvalue. We have hence proved 
the existence of two exact 1/2 eigenvalues for the two 
site problem, cut in half, when the difference between 
the number of negative inversion eigenvalues at 0, 7r is 
±1. 

As in the one-band case, this argument can be ex- 
tended analytically to a system with four sites as shown 
in Appendix [B] indicating that the conclusions hold for 
chains longer than two-sites. 



5. N Occupied bands, Two-Site Problem 

We now show that the two-site problem with n\ neg- 
ative inversion eigenvalues at k — and n 2 negative in- 
version eigenvalues at k = n with a total number N of 
occupied bands contains 2|m — n 2 \ zero modes in the en- 
tanglement spectrum when a real-space cut is made on a 
system with periodic boundary conditions (i.e. there are 
two cuts). The simplest case, which should be obvious 
from our previous examples, is that all the N inversion 
eigenvalues at k = are identical and are the opposite 
of the N eigenvalues at k = tt. In this case, the pro- 
jector at one of the inversion symmetric fc's annihilates 
all the eigenstates at the other inversion symmetric k, 
and the 2N occupied eigenstates of the original two-site 
Hamiltonian are also the eigenstates of the entanglement 
spectrum at fixed eigenvalue 1/2. Due to their orthog- 
onality, they are linearly independent. From here it is 
clear that our formula is correct for this case. 

Now we will prove the more general formula. Let ri\ 
and n 2 be the number of eigenvectors for the —1 eigen- 
value of PqVPo and P^VP^, respectively, and assume 
n\ > n 2 . Recall that K is the number of orbitals per 
site, so V is a KxK matrix acting on C K . V has ±1 
eigenvalues and we denote the invariant subspaces cor- 
responding to the positive/negative eigenvalue by H± 
(H-+H+ = C K ). 

Now the subspaces PqC k and P 1X C K are invariant un- 
der the inversion operation V , and PqC k D H- is pre- 
cisely the subspace spanned by the n\ eigenvectors of 
PqPPq corresponding to its negative eigenvalue. Sim- 
ilarly, P- K C K n H- is precisely the subspace spanned 
by the n 2 eigenvectors of P n VP^ corresponding to its 
negative eigenvalue. Since dim[P C K n H-] = ni and 
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dim[P 7r C A P\H-} = ri2, with n\ > 112, we can always find 
m —n,2 vectors ^ n in PqC k <~)H- that are orthogonal to 
any vector in P^C K PiH- . Since these vectors are in T-L- , 
they are also orthogonal to any vector in P- K C K D H+. 
In other words, \? n 's are nx — n 2 vectors in PqC k per- 
pendicular to all the vectors in P 7T C K . Consequently: 

M„ = i(P + P ff )*„ = itf„, (62) 

for all rii — n-2 vectors Following the same arguments, 
we can find n\ — ri2 vectors in PnC K H H+ that are 
orthogonal to PqC k , and consequently another set of 
m — ii2 eigenvectors with eigenvalue ^. In total, there 
are 2(774 — '^2) robust modes of Cr, at 1/2. For a more 
explicit proof see Appendix \C\ 



6. Extension to higher dimensions 

The extension to higher dimensions is just a matter of 
reinserting the extra momenta which are conserved in the 
presence of the cut. As an example let us consider 2D 
with an entanglement cut parallel to the y-axis so that 
k y is a conserved quantum number. The exact mid-gap 
modes in the entanglement spectrum will exist only at 
inversion symmetric points in the momentum parallel to 
the cut, i.e. k y = 0, it. However, since bands become 
continuous when k y is finely discretized, in the thermo- 
dynamic limit the existence of 1/2 modes at these two 
discrete k-points implies the existence of mid-gap bands. 
The number of bands will be equal to the number of 1/2 
modes. These bands typically disperse away from 1/2 
but do not have to connect with the "bulk" entanglement 
bands at entanglement eigenvalues close to 0, 1 (in the 
special cases of Chern insulators and time-reversal sym- 
metric non-trivial insulators, they do connect and have 
spectral flow, as pointed out previously). However, since 
the modes at 1 /2 are robust upon changes in the original 
Hamiltonian that do not close the band gap, the mid- 
gap bands cannot be entirely pushed to the entanglement 
bulk band edges at or 1. This means the system is a 
non-trivial insulator (i.e. S en t cannot be made to vanish) 
if the number of negative (or positive) inversion eigen- 
values is different between inversion symmetric points. 
We would like to know at which inversion symmetric k y 
the exact 1/2 modes will appear occur. The answer is 
simple: there will be exactly 2\ni — n 2 \ 1/2 modes at 
k y = when the number of negative inversion eigenval- 
ues at (k x , k y ) = (0, 0) differs by \ni — rj-2 1 from the num- 
ber of negative inversion eigenvalues at (k x , k v ) = (n, 0). 
Additionally there will be exactly 2\n' 1 — n' 2 \ 1/2 modes at 
k y = 7r when the number of negative inversion eigenvalues 
at (k x ,k y ) = (0,7r) differs by \n[ — r n! 2 \ from the number 
of negative inversion eigenvalues at (k Xl k y ) = (71", tt). The 
generalization to higher dimensions is the trivial exten- 
sion of this. 



D. Spectral Flow in the Entanglement Spectrum 

As mentioned earlier in this section the two important 
features of the entanglement spectra of inversion symmet- 
ric insulators are protected mid-gap modes, and spectral 
flow. What we mean by spectral flow is a continuous 
connection between the valence and conduction bulk en- 
tanglement bands through the entanglement edge states 
(an example is seen in Fig. [8ji,i ). For a TRI topological 
insulator in 2D and 3D, or for a Chern insulator in 2D, 
the entanglement spectrum mirrors the energy spectrum 
of the open-boundary Hamiltonian. In fact, we have al- 
ready shown an explicit map between the entanglement 
spectrum and the energy spectrum of the open bound- 
ary spectrally flattened HamiltonianPSl This implies that 
if there is spectral flow between the conduction and va- 
lence bands in the energy spectrum then such a flow ex- 
ists in the entanglement spectrum. In fact, this is the 
only case where there is true spectral flow in the en- 
tanglement spectrum. Out of the entire set of inversion 
invariant topological insulators only a small subset have 
spectral flow. Instead most non-trivial systems simply 
exhibit protected mid-gap states (or bands) but these do 
not continuously interpolate between the bulk entangle- 
ment bands. 

There is a nice example which illustrates this di- 
chotomy. Let us consider the 3D strong topological insu- 
lator with both inversion and time-reversal symmetries. 
If we preserve V but break T, spectral flow generically 
disappears from the energy spectrum because gaps are 
opened in the surface state spectrum. The degenera- 
cies which existed in the entanglement spectrum at the 
time-reversal invariant momenta when T is preserved, are 
almost all broken, with the exception of the protected de- 
generacy for states at £ = 1/2. This degeneracy splitting 
breaks the spectral flow in the entanglement spectrum, 
and opens gaps at the TR-invariant momenta as shown 
for a specific model in Fig. [9}i,i. As such, the entan- 
glement spectrum is capable of distinguishing the subtle 
difference between topological insulators with T and V 
symmetry from topological insulators with only V sym- 
metry. 



III. THE LINEAR RESPONSE 

To date, some of the most spectacular features of topo- 
logical insulators are their responses to external fields. 
The Chern insulators exhibit a quantized Hall effect and 
the 3D TRI topological insulators exhibit a topological 
magneto-electric effect. The topological invariants which 
distinguish these states from trivial insulators are directly 
connected with the corresponding response coefficient. In 
fact, a whole ladder of topological responses was uncov- 
ered in Ref. [9j With this precedent one would hope that 
the inversion invariant topological insulators would also 
exhibit some type of defining physical response. How- 
ever, this turns out to be true in only a limited set of the 
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inversion invariant topological insulators. The situation 
is quite varied (remember we only have generic access to 
the information held in the inversion eigenvalues): some 
insulators have unique well defined topological responses, 
some systems can exhibit one of several allowed topolog- 
ical responses, and for others it is unclear if there is any 
topological response at all. We will see examples of all 
three cases in Sec. |IV| In this section though we focus 
on the first case where insulators do exhibit a unique re- 
sponse which is the most interesting physical case. We 
discuss responses in 1,2, and 3 dimensions and then we 
briefly mention how the general pattern might be ex- 
tended to higher dimensions to make contact with Refs. 
|9] and E] in Appendix [K] 



points is the product of the inversion eigenvalues at that 
inversion invariant point: 



(68) 



i=l 



We now turn to the calculation of the polarization 



Pi = — I dkA(k), 

Z7T 



where A(k) is the adiabatic connection: 



A(k) = -i ^2 (ui,k\Vk\ui,k)- 



(69) 



(70) 



A. ID inversion symmetric insulators 

The following discussion applies to a generic one di- 
mensional -RT-band insulator with N occupied bands, and 
with a generic inversion symmetry i.e. 



VU{h)V- x =H(-k), 



(63) 



where the inversion matrix V is unitary and squares to 
the identity: 



V 1 V = 1, VV = l. 



(64) 



As explicitly shown in Appendix [D] the charge polariza- 
tion Pi of a ID insulator behaves as: 



P -> -P 



l +je 



(65) 



under inversion, where j is a gauge-dependent integer. 
This shows that the polarization of ID systems with in- 
version symmetry can take only two values, and e/2, 
modulo a gauge-dependent integer multiple of eP^ We 
will prove that if the xv invariant, defined as the prod- 
uct of all inversion eigenvalues of the occupied bands, 
takes the value 1, then Pi = 0, and if xv — ~ 1> then 
Pi = e/2. More precisely, we will establish that: 



Pi 



2m 



t Log 



iV 



nCi(o)GM 



(66) 



The integer ambiguity of the logarithm is identical to the 
integer ambiguity of the polarization. 

For this, we define the fc-dependent NxN unitary ma- 
trix B(k) 



Bij(k) = (ui- k \V\u 



(67) 



where the indices i and j run only over the occupied 
bands. The inversion eigenvalues Ci(0) and d(ir) coin- 
cide with the eigenvalues of the matrix P(fc), when eval- 
uated at the special inversion fc-points fc=0 and n. It is 
then obvious that the determinant of B(k) at these fc inv 



We will use the following important relation, which is 
proven in Appendix |E} 

A(-jfe) = -A(k) + iTr[P(fc)V fc P t (fc)]. (71) 

The last term can be written in the equivalent form: 

Tr[P(fc)V fc P t (fc)] = -V fe Log[dct[P(fc)]]. (72) 

We can now proceed as follows: 

Pi = i; Jo dk[A(k) + A(-k)] 

= £lfo dk V fc Log[det[£(fc)]], 

with the final answer: 



(73) 



P 



e 
2vri 



Log[det[B(7r)]] -Log[det[B(0)]] . (74) 



This, together with Eq. [68] and the fact that the deter- 
minants can take only the values ±1, so that det[P] 
l/det[P], prove the statement of Eq. 



66 



We mention that similar arguments were used in Ref. [9] 
to classify ID particle-hole symmetric insulators via a Z2 
invariant. In fact, the Zi invariant found there is exactly 
the value of the charge polarization modulo an integer. 
For a ID model with both inversion and particle-hole 
symmetry, such as the ID lattice Dirac model, the in- 
variants coincide. In the non-trivial phase, the ID Dirac 
model exhibits mid-gap energy modes bound to the ends 
of an open chain. The requirement of particle-hole sym- 
metry restricts these modes to lie at zero energy if there 
are an odd number of them. An even number on each 
end is not stable and the degeneracy can be lifted, which 
is another manifestation of the Z 2 nature. The mini- 
mal case is one mode on each end and with particle-hole 
symmetry at half filling one mode is filled and one is 
empty. This leads to an excess charge of +e/2 on the 
side with the filled state and —e/2 on the empty side. 
If we break particle-hole symmetry but keep inversion 
symmetry then both modes can be empty or filled, but 
they are empty or filled together because inversion sym- 
metry connects the two modes. This means that the 
excess charge is either +e/2 on both ends or —e/2 on 
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both ends. However, because of the gauge- variance of the 
polarization this is equivalent to the polarization in the 
particle-hole symmetric case. Thus, both insulators have 
the same topological electric response. One can form a 
complimentary argument by using the effective response 
action for particle-hole symmetric insulators given in Ref. 
H]and[371 



-71 71 



S, 



1 



eff 



dxdtPie^Fau 



(75) 



where F^ v is the field-strength tensor of the externally 
applied electro-magnetic field. The argument for the 
quantization of P\ is as follows. In the partition function 
the phase due to this term is e lS ° ls and under particle- 
hole symmetry Pi — > —Pi. Thus if our system is to be 
particle-hole symmetric we must have e 2jSc// = 1. For 
constant Pi the integral gives 2Tm for integer n and we 
have e 4 ™ lPl — 1 and thus Pi =0,1/2 mod Z in dimen- 
sionless units. Since Pi — > — Pi under inversion symme- 
try the same argument holds and Pi is quantized there 
as well. A similar argument for magneto-electric polar- 
izability of 3D insulators with inversion symmetry was 
given in Ref. [29l As an aside, we recall that in ID 

(2) 

we also had an invariant x-p which helped classify the 
4-band insulator example. We do not know of any re- 
sponse related to this invariant in ID inversion invariant 
insulators. 

Additionally, Eq. [66] can be derived using an alterna- 
tive approach based on a monodromy argument. The 
effect of a magnetic flux $ through a large one dimen- 
sional ring can always be gauged away by a transforma- 
tion: "P— >-e~ l *'I', which is equivalent to a translation in 
fc-space by $. The evolution of the states in response 
to an adiabatically slowly varying magnetic flux can be 
understood from the evolution of the Bloch states in re- 
sponse to an adiabatic translation of fc-spaceP^The mon- 
odromy U(k, ko) describes the evolution of the occupied 
Bloch states when the k of the Bloch Hamiltonian H (k) 
is adiabatically varied, and it is the unique solution of 
the equation: 



^U(k,ko)=i[P k ,d k Pk}U(k,ko), 



(76) 



with the initial condition U(ko,ko)=P ko , assuming that 
we start the evolution from ko. 

The monodromy U(ko, ko) maps the space of occupied 
Bloch states P ka C K at fc into the space of occupied 
Bloch states P k C K at k. Since k = ±ir are the same, 
U(tt, —it) takes the space P- 1V C K into itself, and is a uni- 
tary operator that we call U~, where 7 is one of the paths 
shown in Fig. [3] U 1 gives the change occurring after a full 
quantum of magnetic flux has been pumped through the 
system. If bases {tpi(k)} were pre-chosen for all P k C K 
spaces, the matrix Uij(k) = (ipi(k)\U(k, k )\ipj(ko)} sat- 
isfies the parallel transport equation 




FIG. 3. The adiabatic transport is carried over 71; 



where A(k) is the full non-abelian adiabatic connection 
discussed by Wilczek and Zee in Ref. [35] 

There is a direct relation between the determinant of 
the monodromy and the line integral of the abelian con- 
nection (the trace of the full connection): 



det[U(k f , h)] = exp I / Ti[A(k)}dkj . (78) 

Indeed, working with predefined basis sets for P k C K and 
breaking the interval kf, fcj in small subintervals kf, k n , 
. . ., hi, we have, up to second order corrections: 



U{k f ,ki) = U(kf,k n )U(k n ,k n -i).-.U{ki,ki) 



(k f - k n )A(k n )) 



(ki-ki)A(ki)). 



(79) 



Taking the determinant on both sides and using some 
elementary identities, we obtain 

det[U(k f , h)] = (1 + (k f - fc„)Tr[i(fc„)]) 

. . . (1 + (h - ki)Tr[A(ki)]) (80) 
= e (fe / -fe»)Tr[4(fe n )] _ _ _ e (k 1 -fc < )Tr[A(fc i )] j 

from which the identity of Eq. [78] follows. The identity 
is valid in higher dimensions too. 

Assuming k =—7r, a conjugation of Eq. [76] with the 
inversion operator V gives: 

d k {TU(k, -^V- 1 } = i[P- k ,d k P^ k ]VU(k, -n)V- 1 , 

with the initial condition VU(-ir, -tt)V^ 1 =P(tt). This 
is just the equation for U(— k, it), which shows that 
VI] (k, —TiyP~ x coincides with U(—k, tt). Equivalently we 
can think that V sends 7 into 7' in Fig. [3] Now obviously 
U 1 U 1 ' equals the identity, therefore: 



det^PKyP- 1 ] = 1. 



(81) 



Using the elementary properties of the determinant, we 
conclude that det[C/ 7 ] 2 =l, hence det[C/ 7 ] can take only 
two values: 



mU(k) 



A(k)U(k), 



(77) 



det[t/ 7 ] = ±1. 



(82) 
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The following calculation will show that the cases 
det[£7 7 ]=±l correspond to Pi=0 and Pi = § (mod Z), re- 
spectively. Indeed: 



det[P 7 ] = det[J7(7r, 0)U(0, -tt)] 
= det[C/( 7 r,0)7'C/(0,7r)7'- 1 ] 
det[^(7r, 0)P PP P(0, ttJP^ -1 ^ 



(83) 



In the last line we have inserted the projectors Po j7r to see 
explicitly the spaces on which P is acting. Using again 
the elementary properties of the determinant and the fact 
that U(tt, 0)U(0, 7t)=1, we obtain: 



det[[> 7 ] = det[P PP ] detlArP" 1 ^] 



N 



det[C/ 7 ] = TT Ci(0)CiW- 



and thus 



Pi 



2™ 



Log det[£7 7 



(84) 



(85) 



(86) 



(cf. Eq. 78 1 



For a topological insulator, no matter what definition 
one uses, it is always the case that, when the hopping 
terms between the neighbors are adiabatically turned off, 
that is, when one takes the atomic limit, the insulating 
gap of the system closes at some point in the process. One 
can investigate this issue using the inversion eigenvalues 
directly, as we have focused on, but here we see how the 
physical response enters into the picture. Note that the 
inversion eigenvalues can be easily computed for simple 
models, but may not always be available, especially for 
complex materials. Pi or U 1 are physically measurable, 
so they can provide physical signatures of the non-trivial 
state. In the atomic limit, the bands have no disper- 
sion, so in this limit C/ 7 is just the identity matrix. If 
Det[C/ 7 ]=— 1, it is obvious that U 1 cannot be smoothly 
connected to the identity and the insulator is topological. 
However, we also must note that if Dct[P 7 ]=l, it is not 
necessary that the insulator is trivial. 

We can refine the investigation by asking when can £7 7 
be smoothly connected to the identity? For this, let us 
look again at Eq. [83] 



P 7 = U{ir,0)PoVPoU(p,Tr)P*PP v . 



(87) 



The first term, U(ir, 0)P o VP o U(0, tt), is just P PP par- 
allel transported from k=0 to k=ir. The parallel trans- 
port does not alter the eigenvalues of P(0)PP(0), which 
are pinned at ±1 (recall V 2 —l). The eigenvalues of 
P(tt)VP(tt) are also pinned at ±1. So what we have 
in Eq. [87] is a product of two operators with eigenvalues 
pinned at ±1 and because of that the eigenvalues cannot 
change under any smooth deformation of the Hamilto- 
nian that keeps the insulating gap open and preserves 



inversion symmetry. Now, if U 1 can be deformed into 
the identity, then U(ir, 0)P o PP o P(0, n) can be turned 
into the inverse of P^VP^ and this requires that PoPPo 
and P^VP-k have identical eigenvalues, counting the de- 
generacy too. The conclusion is: if the set of inversion 
eigenvalues of 



P PP and P W PP W 



are not identical, then U 1 cannot be connected to the 
identity. Since the eigenvalues are restricted to just ±1 
values, then the following integer: 



X : 



Tr PqPPq - P.PP 



■A 



(89) 



tells how many eigenvalues are different for the two ma- 



trices in Eq. 88 To reach the atomic limit, we need to flip 



X inversion eigenvalues, and that will require a minimum 
of x S a P closings. Unfortunately, we were not able to 
find an expression of the topological invariant x solely in 
terms of the monodromy P 7 , but we know there are pre- 
cisely x topological obstructions when trying to connect 
the monodromy to the identity. The topological invari- 
ant x also gives the number of robust edge modes seen 
in the entanglement spectrum on a single cut as shown 
in Sec. [H] 



B. 2D inversion symmetric insulators 

The physical response of 2D inversion symmetric in- 
sulators is much richer than that of ID. Based solely on 
the inversion eigenvalues, one can define several invari- 
ants, the first of which is the isotropic extension of xp to 
2D i.e. 



Xv 



(90) 



where fc inv runs over all four inversion invariant fc-points. 
We show that 



Xv = (-1) 



Ci 



(91) 



where G\ is the first Chern number of the ground state. 
Thus if xv is negative the system must be in a quantum 
Hall state, and if it is positive it is in a state with an 
even Chern number which can be zero. Thus, only if it is 
negative are we sure it is in a topological insulator state. 

We will prove the statement of Eq.[9l]m two ways, first 
using a band crossing argument and then a monodromy 
argument. Let us begin by assuming we are in a trivial 
insulator state in an atomic limit with N occupied bands 
and that we have inversion symmetry. We can reach 
any non-trivial topological insulator state from this limit 
through a series of Hamiltonian deformations that will 
lead us through band crossings. Our assumption of an 
atomic limit implies that xv = +1 initially. If we want 
to generate xv = — 1 we need to have an odd number 
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FIG. 4. The two paths 7 and 7 used in the monodromy 
argument in 2D. 



of band crossings between bands with opposite inversion 
eigenvalues. Assume we have an odd number of such 
crossings. This implies that there must be an odd number 
of crossings at the inversion invariant points. This is true 
since crossings occurring at non-invariant k are accompa- 
nied by a crossing at — k always giving an even number of 
eigenvalue switches which will not affect the value of \v ■ 
Thus we only need to consider the odd number of cross- 
ings occurring at the invariant momenta. The generic 
Hamiltonian of each crossing between opposite inversion 
eigenvalue states near an inversion invariant momentum 
is H = p\o x + P2O2 + m e ffa 3 where {pi,P2) is the mo- 
mentum away from the inversion invariant point, m e ff is 
a term parameterizing the distance to the band crossing, 
and cr 3 is the inversion operator projected onto the two 
crossing bands. Exactly at the inversion invariant mo- 
mentum the Hamiltonian reduces to m e ffa 3 as it must 
in order to commute with V there. As the crossing occurs 
m e f f switches sign and an inversion eigenvalue of the oc- 
cupied band is changed. We note that this Hamiltonian 
is a 2D massive Dirac Hamiltonian which is switching the 
sign of the mass. At such a crossing, the Chern number 
changes by ±1 and thus we see that going through an 
odd number of crossings changes the parity of the Chern 
number. Thus if \v = — 1 the parity of the Chern num- 
ber is odd since the parity is even in the initial atomic 
limit i.e. C% — 0. The inversion eigenvalues thus give us 
a rough way to characterize the quantum Hall effect in 
an insulator. 

The corresponding monodromy argument proceeds as 
follows. We consider the monodromy corresponding to 
the path 01234560 in Fig. |1J starting from the middle 
of the Brillouin zone and continuing on its rim. This 
path can be also be viewed as the composition 7+7' of 
the paths 7=012340 and 7'=045610, and the monodromy 
corresponding to 7+7' can be written as a product of 
partial monodromies: 



U, 



7+7' 



U iU 16 U 65 U 54 U 43 U 32 U 21 U 10 . 



(92) 



Since the products U10U01, C/16^54, ^65^32 and U43U21 
are all equal to the identity, taking the determinant of 



Eq. 92 leads us to conclude that det[{7 7+7 ']=l. This is 



e 2iu-Ci ^ see gq_ 78 1 an d the first Chern number is an 
integer. 

The next observation is that inversion sends 7 into 7' 
and consequently: 



det[U 1+Y ] = detp^VUjV] = det[U 7 } 2 



(93) 



The conclusion is that the determinant of U 1 can only 
take the values ±1. Since the path 7 encircles half of the 
Brillouin zone, and the adiabatic curvature is symmetric 
when k — > —k, it is also true that det[C/ 7 ] = e J7rCl (see 



Eq. 78) 



We now take a closer look at U^,. Since inversion sends 



10 into 40 and 20' into 30', we have: 



Uj — C/o4^ r 43^30'^0'2^ r 2lt^lO 

= VUoiVU 12 VU20'VUa'2U2iU w 



(94) 



Inserting the appropriate projectors to specify explicitly 
on which spaces are the V operators acting, taking the 
determinant and using its elementary properties together 
with the fact that UijUji equals the identity, we obtain: 

(-l) c i = det[l> 7 ] = det[P , PPo,o] det[P , w pP 0j7r } 
x [det{P n>v VP^ n ] det[P w , PP T ,o], 

(95) 

which completes monodromy argument for Eq. |91| 

We will briefly mention two other interesting inversion 
invariants one can define in 2D. The first is an anisotropic 
invariant defined by taking the product of the inver- 
sion eigenvalues at only two invariant momenta in the 



not surprising and is related to the fact that det[C/. 



7+7' J 



2D Brillouin zone. We provide an example in Sec. IV 
with XT = +1 but where this anisotropic invariant is 
non-trivial. This invariant has interesting implications 
for the charge polarization, but there are some subtleties 
that we will illustrate. To begin, assume we have an in- 
version invariant, insulator Hamiltonian H(k x ,k y ) with 

N occupied bands such that YliLi C*(0> 0)Ci(0; = — 1- 
This implies that the ID Hamiltonian H(Q,k y ) has a 
charge polarization Pi = e/2 since this restricted Hamil- 
tonian is inversion invariant. Now we want to know if 
this non-trivial anisotropic invariant is enough to specify 
the full polarization of the entire 2D system. The an- 
swer is no. To see why, we slightly deform P(0, k y ) away 
from the k y axis. The ID Hamiltonian H(Sk x ,k y ) is not 
generically invariant under inversion symmetry because 
it gets mapped onto the Hamiltonian H(—Sk x ,k y ) and 
thus the polarization does not have to remain quantized 
with value e/2p3ln fact, by the time we have deformed 
all the way to the Hamiltonian H(jr, k y ) the polariza- 
tion, which must again be quantized since this Hamilto- 
nian does have inversion, can be completely different. So 
while we can think of H{k x , k y ) as a gapped interpolation 
between H(0, k y ) and H(ir, k y ), inversion symmetry, and 
thus the value of the polarization, is not preserved along 
the interpolation. Intuitively this makes sense because 
it is exactly when the polarization changes its quantized 
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value that the system has an odd Chern number. An odd 
Chern number is allowed because a quantum Hall effect is 
not forbidden by the requirement of inversion symmetry 
Thus, if we have 2D inversion symmetry the anisotropic 
invariant does not determine the 2D charge polarization. 

As an aside, since we clearly know why inversion sym- 
metry fails, we immediately know how to fix the problem. 
We fix it by requiring a reflection symmetry (i.e. parity 
symmetry) about an axis, instead of inversion. With- 
out loss of generality we can have a reflection symmetry 
M such that MH(k x ,k y )M^ 1 = H(k x ,-k y ). We will 
see that the reflection eigenvalues will specify the po- 
larization in the y-direction (for this choice of reflection 
symmetry). First, we see that the inversion invariant 
momenta are also reflection invariant momenta. Thus 
[H(ki nv ),M] = 0, and we can label the states at these 
points by their reflection eigenvalues which we will also 
call Q(k inv ). Now let us assume the same setup as the 

previous paragraph with Yli=i Ci (0, 0)C,(0, tt) = — 1. The 
point is that now when we adiabatically deform away 
from the k y axis the ID Hamiltonian H(5k x ,k y ) is in- 
variant under reflection. Additionally, the y-component 
of the polarization is quantized as long as reflection is 
a good symmetry. Thus, H(k x ,k y ) is a gapped inter- 
polation along which reflection symmetry is preserved 
and thus the y-component of the polarization is fixed 
and quantized to e/2 for each ID Hamiltonian. This ar- 
gument immediately implies that J\i=i di^: 0)0 ( 7r i 7r ) — 
— 1. Thus the parity of the Chern number is always even 
when reflection symmetry is required. In fact, it always 
vanishes because the quantum Hall effect is incompat- 
ible with reflection symmetry. Notice that the reflec- 
tion eigenvalues do not uniquely specify the polariza- 
tion in the x-direction since the two eigenvalue condi- 
tions can be satisfied by choosing Yl i=1 C»(0, 0)Q(ir, 0) = 

rii=i Ci(0) 7r )Ci( 7r : 7r ) = il- Thus, the anisotropic invari- 
ants in the x-direction can take either value. 

The other invariant one can define is the isotropic ex- 

(2) 

tension of x-b to 2D. If \v = +1 and the product over 
the inversion eigenvalues at every invariant momentum is 

(2) 

separately trivial then \ v is well-defined because there 
are an even number of negative inversion eigenvalues at 

(2) 

each invariant momentum and by definition x v IS the 
product over half of those negative eigenvalues. This is 
the Z2 topological invariant which indicates a quantum 
spin Hall effeclplwhen an insulator has both inversion and 
time-reversal symmetry. Unfortunately this topological 
invariant does not yield a unique topological response. 
We can understand this in several ways. By explicit con- 
struction take two decoupled copies of the QAHE state 
each with C\ = 1 to give an IQHE with C\ — 2 or 
take two decoupled copies of QAHE one with C\ = 1 
and the other with C\ = —1. The first system breaks 
time-reversal and still gives an IQHE while the second 

preserves time reversal and will not give a quantum Hall 

(2) 

effect. These systems both have xv — +1 but x-b — — 1- 
We can immediately see why this invariant does yield a 



unique response in the presence of an additional time- 
reversal symmetry because this requires the total Chern 
number to vanish which only leaves the possibility of a 
quantum spin Hall state. Thus it is time-reversal symme- 
try which restricts the allowed physical response. We can 
understand this by a simple symmetry argument as well. 
The quantum spin Hall effect is even under both inversion 
symmetry (x,y) — ► (—x,—y) and parity (x,y) — > (x,—y) 
We see that in even space dimensions inversion symmetry 
is just a rotation. The quantum Hall effect is even under 
inversion but odd under parity. Thus having a quantum 
Hall effect is compatible with inversion symmetry but not 
parity (and not time-reversal) . This is why the inversion 
classification does not distinguish between the doubled 
quantum Hall state and a quantum spin Hall state. This 
type of ambiguity exists in every dimension where one 
can define a Chern number invariant. However, as we 
saw with the polarization, if we consider the eigenvalues 
of a reflection symmetric Hamiltonian we can eliminate 
the possibility of a non-zero Chern number thus leaving 
a quantum spin Hall state as the alternative. 

C. 3D Inversion Symmetric Insulators 

We show that the isotropic extension of xv in 3D 

XV = 1] Ci(*tor), (96) 
/ci nv ;zGocc. 

where the product runs over all 8 inversion symmetric 
fc-points, is always trivial. We can first see this by using 
a band crossing argument. Start from the atomic limit 
in which bands have identical inversion eigenvalues at 
all inversion symmetric points. For xv to be non-trivial 
(equal to —1) there has to be an odd number of band 
crossings between bands of opposite inversion eigenval- 
ues at the inversion invariant momenta when beginning 
from this trivial limit. Without loss of generality let us 
consider one crossing between two bands with opposite 
inversion eigenvalues at all inversion symmetric points in 
the Brillouin zone. The crossing can happen at either 
an inversion symmetric point (by tuning one parameter) 
or at a non-inversion symmetric point k in the Brillouin 
zone. In the latter case, there are actually two cross- 
ings at k and —k because of inversion symmetry. To 
initially close the gap between the two bands one needs 
a quadratic touching in at least one of the directions, 
otherwise we would be creating a nonzero Chern num- 
ber Fermi surface (after gap closing) out of a zero Chern 
number surface(before the closing, in the trivial limit). If 
the crossing starts at an inversion symmetric point, the 
quadratic touching and the gap reopening at the inver- 
sion symmetric point will switch the inversion eigenvalues 
of the bands at that point and make xv = — 1. However, 
the system will no longer be an insulator: it will have 
two crossings somewhere else in the Brillouin zone. That 
is, although the gap has opened at the inversion sym- 
metric point, the gapless points have been moved away. 



18 



The quadratic touching effectively splits into multiple 3D 
Dirac points. In the generic situation, there are two gap- 
less points in the Brillouin zone, with relative position 
fixed by inversion, and in order to gap the system, we 
need to annihilate the two Dirac points. Note that a 3D 
Dirac point is locally stable even if inversion is not pre- 
served. On a 2D surface surrounding a 3D Dirac point 
(Hi oca i(k) = kiAijVj, (i.j = 1,2,3)), the Chern number 
is C — sgn(det(A,-j)). The degeneracy points are stable 
unless two Dirac points of opposite Chern numbers anni- 
hilate. Inversion symmetry forces the points at k, — k to 
have opposite Chern numbers. So, by inversion, annihila- 
tion of only two Dirac points in the whole Brillouin zone 
can only happen at another inversion symmetric point, 
in which case inversion eigenvalues are switched again to 
give xv — 1 f° r an insulator. If the gap first closes at 
a non-symmetric point k, we find the same end result: 
generically two Dirac points are created close to k and 
two close to —k. They can annihilate in pairs, always 
switching an even (possibly zero in this case, since the 4 
dirac points can annihilate two by two at non-inversion 
symmetric points in the BZ) number of inversion eigen- 
values. 

Another similar way of understanding the 3D band 
crossings is the following: since we are considering a 
two-band crossing in 3D there are three varying param- 
eters and we cannot find a gapped phase with xv = — 1, 
though it is possible to find a gapless phase. If we 
have an even number of band crossings then we can 
open a gap but this means that two inversion eigen- 
values are switched leaving xv = +1- There is also a 
deeper reason why this cannot be done. Let us look at 
the simplest case of a single two-band crossing at the F- 
point. The effective Hamiltonian can be put in the form 
H e ff = piG 1 + p2<7 2 + p 3 cr 3 which is a (chiral) Weyl- 
fermion Hamiltonian. From the Nielsen-Ninomiya no-go 
theorem this type of Hamiltonian cannot arise without 
a partner fermion with the opposite chirality.^l Thus all 
two-band crossings must generically occur in pairs and 
there cannot be an odd number of negative parity eigen- 
values in an insulating system. Conversely, if no states 
at the inversion invariant momenta cross the Fermi-level 
and xv — ~ 1 we immediately know that the system con- 
tains a gapless point(s) somewhere in the Brillouin zone. 

We now provide an alternate proof that the product 
of all the inversion eigenvalues of all the occupied bands 
must be positive in 3D. This will prove useful when con- 
sidering possible 3D QHE states. Assume a gapped in- 
sulator Hamiltonian H(k x , k y , k z ) which has N occupied 
bands and is invariant under inversion symmetry. We 
can take the plane k z — 0, and think of H(k x , k y , 0) as a 
2D inversion symmetric Hamiltonian. We have already 
proved that the inversion eigenvalues of a 2D inversion 
symmetric Hamiltonian are related to the Chern number 
i. e.: 

nti Ci(0, 0, 0)Ci(7r, 0, 0)0(0, tt, 0)0 (tt, tt, 0) 



The same thing is true for the inversion symmetric 2D 
Hamiltonian H(k x , k y ,ir): 

UiLl 0(0, 0, w)Q(lT, 0, 7T)Ci (0, TT, 7r)C 4 (7T, TT, 7r) 

= (_i) c 'i(^( fc »> fc v> ,r )), (98) 

Hence 

11=1 Ci(0, 0, 0)0(71, 0, 0)0(0, tt, 0)0 (tt, tt, 0) 

• 0(0, 0, 7r)0(7T, 0, 7r)0(0, TT, 7r)0(7T, 7T, 7t) 

= (_l) c i( j& (**.*v.0))+Ci(^(k«,fe v ,ir)) < ( 99 ) 

Since H(k x ,k y ,k z ) is gapped due to our assumption of 
an insulator, we can think of it as an adiabatic interpo- 
lation between H(k x ,k y ,0) and H(k x ,k y ,Tr) by varying 
the parameter k z . Since this interpolation preserves the 
U(l) charge conservation symmetry i.e. there is no su- 
perconductivity, it means that the Chern number cannot 
change from k z — to k z = tt. Thus, 

Ci{H{k x ,k y ,Q)) = Ci{H{k x ,k y ,TT)) = Ci. (100) 

Hence 

J] 0(fc,nv) = (-l) 2Cl =l. (101) 

ki nv ;i£occ. 

1. Anisotropic Invariants and the 3D Quantum Hall Effect 

We have seen that the isotropic invariant in 3D, which 
is constructed by multiplying the inversion eigenvalues at 
all invariant momenta, is always trivial for an insulator. 
However we can form anisotropic invariants by consider- 
ing the products of inversion eigenvalues over planes or 
lines of the Brillouin zone which are mapped onto them- 
selves under inversion. 

First consider a plane in the 3D Brillouin zone which is 
mapped onto itself under inversion. To be explicit, take 
the plane k z = tt. If the product of inversion eigenvalues 
of all the occupied bands in that plane is 

N 

[]0(0,0,7r)0(7r,0,7r)0(0,7r,7r)0(7r,7r,7T) = -1 (102) 
i=i 

we have a 3D QHlP^ with 3D Hall conductance 

a xy = odd integer x — (103) 

where c is lattice constant in the z direction. From our 
above proof we know that the product of eigenvalues in 
the k z — plane must also be —1. The proof of the 
3D QHE is simple: H(k x ,k y ,k z ) can be thought of as 
an adiabatic continuation of H{k x ,k y ,Tr) (since we as- 
sumed it to be an insulator) As such, each k z plane 
H(k x ,k y ,k z = constant) has an odd integer QHE. Mul- 
tiplying by the momentum, we get the above 3D Hall 
conductance. It is important to note that this argument 
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depends crucially on the fact that the Chern number re- 
mains unchanged when performing adiabatic deforma- 
tions as long as charge conservation symmetry is pre- 
served (i.e. wc don't allow superconducting perturba- 
tions) . 

In general the 3D quantum Hall effect directions can 
be inferred from the eigenvalue formulas. The general 
formula for the 3D Hall conductance in terms of the in- 
version eigenvalues is 



c Q/ 3_l 7 = G 7 (2n + 1/2 

-lu n 



Q{kinv))- (104) 



zGocc. ki nv Gplanc_LG-y 

This expression above gives the 3D Hall conductance as a 
product of the inversion eigenvalues in a plane a/3 perpen- 
dicular to the y(= x, y, z) direction. G 7 is the reciprocal 
lattice vector in the 7 direction and we have left out the 
units of e 2 /h. If the product over all the bands of the 
inversion symmetric eigenvalues in the a/3 plane is —1 
we can see that the Hall conductance is an odd integer 
in that plane and hence cannot be zero. If the product 
is +1 then the 3D Hall conductance is even and can be 
zero. 

Finally, just as in 2D we can consider the product of 
inversion eigenvalues along a single inversion invariant 
line in the 3D Brillouin zone. Again inversion symmetry 
does not allow us to determine the polarization but a 
reflection symmetry about a plane allows us to specify 
the polarization perpendicular to that reflection plane. 
The argument is basically the same as in 2D so we did 
not include it here. The result, however, is explicitly 
illustrated with the 3D dimer model shown in Sec. IIVI 



2. Topological Metal State 

As a corollary to the above, when the product over the 
inversion eigenvalues of several bands at all points in the 
Brillouin zone is negative: 



Xv 



ki nv ;i£occ. 



the system is a metal protected from opening a gap from 
infinitesimal perturbations. We can gain some intuition 
about why this metal state exists by looking at the ef- 
fective Hamiltonian near a band crossing between states 
with opposite inversion eigenvalues at an inversion invari- 
ant momentum. In the most generic case, in which we do 
not have any extra point-group symmetries the effective 
Bloch Hamiltonian expanded for small k near ki nv : 



H(k) = (M + A lj k l k. ] )cr z + kiB ia a a 



(105) 



where i,j — 1,2,3 while a = 1,2. Notice that since the 
bands have opposite parity, the mixing elements between 
them have to be odd in k. The Hamiltonian reduces to 



Ma z at ki nv and thus the sign of M dictates the occu- 
pied inversion eigenvalue. When M is tuned thru zero, 
we have a phase transition with an eigenvalue switch be- 
tween + and — . For a given A^ matrix, notice that on at 
least one side of the switch, we must have a gapless phase: 
if M > and Det(Aij) < we have a metallic phase, or 
if M < and Det(Aij) > we have gapless Dirac points 
away from ki nv . There are of course two Dirac points, 
which can then annihilate at another ki nv point with ± 
eigenvalues and can re-open the gap and give \v = +1- 
Because of the different eigenvalues under inversion, we 
only need to tune one parameter to get a band crossing 
in this case. 



3. Magneto-electric Polarization and Inversion Symmetry 

Although the first 3D isotropic invariant we mentioned 
is always trivial, Ref. [29] argues that inversion invari- 
ant insulators in 3D, which come from strong topologi- 
cal insulators with softly broken time-reversal invariance, 
can support an isotropic topological magneto-electric re- 
sponse (i.e. a 6 = it vacuum). Their argument uses the 
transformation properties of the effective response action 
which we recount here. The topological response action 
in 3D for an insulator coupled to an electromagnetic field 
is 



Seff [A,] = 



2nh 



d 4 x9E ■ B 



(106) 



where E, B are external applied electric and magnetic 
fields, and 9 is an intrinsic quantity which is proportional 
to the magneto-electric polarizabilityP For translation- 
ally invariant systems the magneto-electric polarizability 
for time-reversal invariant insulators is 



Pi 



1 



2tt 32tt 3 



d 3 ke ijk Tr 



AsFi, 



2i A A A 



(107) 

where Ai(k) is the non-abelian adiabatic connection, and 
Fij(k) is the non-abelian adiabatic curvature. Under- 
time-reversal symmetry P3 — > — P3 and thus, for time- 
reversal invariant insulators P3 = 0, 1/2 (in units of 2tt). 
P3 is not a gauge invariant quantity and only defined 
modulo an integerP Note that under time-reversal B 
does not change in the effective action since it is an 
externally applied field. Only intrinsic quantities such 
as P3 get acted upon with time-reversal. Using these 
two values of P3 we can physically define a non-trivial 
time-reversal invariant topological insulator as one with 
P3 = 1/2. It has been shown that this definition is equiva- 
lent to th e band structu re definition of strong topological 
insulators! 7 * 9 * 13 * 14 * 42 * 43 ^ In the presence of time-reversal 
and inversion symmetries there is an elegant topological 
invariant one can define 



■1) 



2F 3 



n n d{ka) 



(108) 
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which is the product of the inversion eigenvalues at every 
invariant momenta in the 3D Brillouin zone for half the 
occupied bands. Since we have time-reversal, half the 
bands just means one band out of each Kramers' pair. We 
provide a new physical proof of this equation in Appendix 

El 

The additional insight of Ref. HH1 is that P3 is also odd 
under inversion symmetry. This means that the values 
of P3 are still quantized to be 0, 1/2 even when the sys- 
tem does not have time-reversal symmetry but only in- 
version. We find that this argument holds for topological 
responses in all even spacetime dimensions (see Appendix 
|k| . Thus, there are inversion symmetric topological insu- 
lators in 3D where P3 = 1/2, analogous to the case in II? 
where Pi = e/2. This argument is an indicator that inver- 
sion symmetric insulators can support non-trivial topo- 
logical states with interesting response properties. Now 
we can ask the question, is Eq. |108| still valid when only 
inversion symmetry is preserved? If we only have inver- 
sion symmetry and no time-reversal symmetry, then this 
formula continues to apply if we can adiabatically con- 
nect the system to the T and V invariant limit without 
breaking inversion symmetry. However this is not the 
only case, and we will show exactly when the inversion 
eigenvalues in 3D indicate a non-trivial magneto-electric 
polarizability protected by inversion symmetry in the fol- 
lowing section. Our arguments use many of the results 
discovered in the previous sections. 



4- Magneto-electric Polarizability for Inversion Invariant 
Insulators 



We begin by reintroducing the unitary matrix 



By (ft) 



V\u 



(109) 



where it^fc is a Bloch function with i, j labeling which 
occupied band and k is the Bloch momentum. An im- 
portant property of the matrix B is 



B(-k) = B f (fc) 



(110) 



which is true because V is unitary and squares to the 
identity matrix. This means that at the inversion sym- 
metric points, the matrix is real and symmetric. For 
inversion symmetric insulators, we prove in Appendix [F] 
that the non-Abelian adiabatic connection satisfies: 



ii(-fe) = -Bi i (fc)B t +iB(k)V l B\k) 



(III) 



which implies that the adiabatic curvature is gauge co- 
variant via: 



F ij (-k)=B{k)F ij (k)B^(k). 



(112) 



From this the magneto-electric polarizability is easy, but 
tedious (see Appendix [f]) , to obtain: 

1 



2P3 



24tt 2 



d 3 ke ijk Tr[(B(k)diB\k)) 



This proves that P3 is either integer or half-integer de- 
pending on whether the RHS (which is an integer winding 
number) is even or odd. Since P 3 is itself defined mod 1 
it means that P3 defines a Z2 classification that indicates 
a non-trivial topological response. 

With time-reversal symmetry, states pair up in 
Kramers' doublets and are degenerate at time-reversal 
invariant points (same as inversion symmetric points). If 
we break all accidental degeneracies in a time-reversal in- 
variant system we are still left with an even number of 
occupied bands - they are degenerate in pairs at time- 
reversal symmetric points. The B{k) matrix then, at 
most, decomposes into diagonal blocks of decoupled U(2) 
matrices. The winding number can be non-trivial in this 
case because the homotopy group n 3 (U(2)) = Z (the 
Z2 nature appears because P3 is defined as the winding 
number mod 2). Once we lose time-reversal symmetry 
the bands do not have to be degenerate at the ki nv . In 
fact, once time reversal is softly broken, the bands at ki nv 
experience eigenvalue repulsion: before TR was broken, 
the inversion eigenvalues of the Kramers' doublets had 
to be identical. As such, it seems that once time-reversal 
is broken we can completely separate the occupied bands 
and isolate them from each other at all points in the Bril- 
louin zone. This would imply that the matrix B could be 
reduced to an N x N diagonal matrix with U(l) phases 
on the diagonal. In other words, B(k) will map from 3D 
momentum space into U(1)® N . Since 7^(^(1)®*) = 
this implies that the winding number in Eq. |113| always 
vanishes (note that by making this statement we are im- 
plicitly assuming that all of the mappings are smooth 
which we will come back to later) . This trivial reasoning 
would seem to imply that we cannot get a topological 
insulator with inversion symmetry. Fortunately, this line 
of reasoning fails because of the non-trivial global con- 
straint that the product of all the inversion eigenvalues 
must be positive. Before we deal with the effects of the 
constraint we draw several conclusions about a Hamil- 
tonian with only one occupied band. The non-Abelian 
form of the winding number implies that with only one 
occupied band P3 = 0. Thus, no matter what the inver- 
sion eigenvalue content, there is no contribution to P3 
from a system with a single occupied band. This also 
holds true for a band that can be completely separated 
and untangled from all the other bands at all points in 
the Brillouin zone. Such isolated bands do not contribute 
to a non-trivial P 3 . This is in contrast to the statement at 
the end of Ref. HUwhich seems to state that a non-trivial 
P3 only requires an odd number of pairs of negative inver- 
sion eigenvalues. As a counter-example, a single occupied 
band can have a single pair of negative parity eigenval- 
ues (which occur at different ki nv ) and it has vanishing 
P3. We give such a Hamiltonian in Eq. 138 We will see 



■{B{k)d 3 B\k)){B(k)d k B\k))\. (113) 



that the important thing to consider is pairs of negative 
parity eigenvalues at the same ki nv . 

The global constraint on the inversion eigenvalues is 
crucial for our discussion. We first give an explicit ex- 
ample to gain intuition about its importance. Let us 
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assume that N = 2 and that the occupied bands have 
Ci(0, 0, 0) = C2(0, 0, 0) = —1 and all other inversion eigen- 
values at all other k inv are +1. The naive reasoning from 
above implies that by perturbing this Hamiltonian while 
preserving inversion symmetry we should be able to sep- 
arate these bands so that they are isolated with no in- 
termingling degeneracies. However, if we could make the 
two bands non-degenerate everywhere we would contra- 
dict the constraint xv = +1- This is easy to see because 
if we could separate the two bands we could consider 
a different insulating ground state with only one of the 
previous two bands occupied and then the product of the 
inversion eigenvalues of the single occupied band would 
be negative. But we have proved that this is not possi- 
ble, so we cannot make the bands fully non-degenerate 
over the entire Brillouin zone. Alternatively, if we could 
separate the bands, we would have a single band that has 
eigenvalues — h ++ in the k z = plane and + + ++ on 
the k z — 7r plane. If we consider the 3D Hamiltonian as 
a gapped interpolation between these two planes then we 
have adiabatically connected 2D Hamiltonians with odd 
and even Chern number respectively. This cannot hap- 
pen and is thus another contradiction. Hence, we can be 
sure that the two occupied bands are degenerate at at 
least two points in the Brillouin zone (by inversion sym- 
metry). These two points are exactly enough to transfer 
an even Chern number from one plane to the other which 
fixes the Chern number issue. Due to this topological ob- 
struction, the B(k) map cannot be reduced to a map from 
T 3 into E7(1)®I7(1), but instead it is a true map from T 3 
into U (2) and therefore its winding can take non-trivial 
values. Additionally, the winding number for each U (2) 
block is additive due to the trace in the winding number. 

The conclusion at this step is that, in general, the B(k) 
matrix can be decomposed into one and two dimensional 
diagonal blocks. The one dimensional blocks have no con- 
tribution to the winding number. The two dimensional 
blocks come from paired bands that cannot be disentan- 
gled. Following the previous arguments, one can see that 
this happens only for the following inversion eigenvalues 
patterns at the eight fc inv (here we restrict the projector 
to just these bands): 



-+ + +++ h++- 

-+ + +++' h++- 



-+- 
-+- 



Note that the pattern (_) at the same fc inv is excluded 
in all these cases. The reason is because, if the pattern 
(i) shows up at a fc inv , then it will necessarily show up 
at another fc inv in order to make xv—^- In this case, a 
simple exercise will show that the inversion eigenvalues 
can always be separated in two groups with the product 
of inversion eigenvalues in each group equal to one. Thus, 
there is no topological obstruction and the bands can be 
untangled. The same thing happens for all the other 
possible inversion eigenvalues that are not listed above. 
We will discuss more about this below and we now return 
to the calculation of the winding number. 

We first discuss the winding number in the continuum 
(sphere) and then focus on the lattice (torus). We know 
from Ref. [TT] that when considering the isotropic topo- 



logical invariants such as P3 it is sufficient to think of 
momentum space as a sphere S 3 instead of the Brillouin 
zone torus T 3 . Allowing for the full torus structure gives a 
rich set of anisotropic states which we will consider later, 
but for now we assume that the momentum space topol- 
ogy is spherical. This effectively reduces the number of 
invariant momenta we need to consider to just two: the 
origin and the 'point at infinity' Equivalently we could 
think of the torus with unrestricted inversion eigenval- 
ues at ki nv = (0, 0, 0) but with the inversion eigenval- 
ues at all the other k inv constrained to be equal band 
by band. Now consider a Hamiltonian with N occupied 
bands (again N does not have to be even, a crucial dif- 
ference with the time-reversal case). Since we are in 3D 
we know that the product of all the inversion eigenval- 
ues must be +1. This means that there can only be an 
even number of inversion eigenvalues that are different 
between k = and k — 00. For example, the number of 
negative eigenvalues at k = 00 must have the same parity 
(i.e. even or odd) as the number of negative eigenvalues 
at k = 0. It is clear that, by exchanging parity eigenval- 
ues between bands, either at k = 0, 00, we can split the 
bands into two classes: (i) bands with \p = +1 an d (ii) 
pairs of bands with xp — ~ 1 f° r each band. We see case 

(i) when the eigenvalues match at k — and 00 and case 

(ii) when they are opposite. For case (i) the bands can 
be isolated from each other, but in case (ii) they must 
generically be in tangled pairs where the xp = +1 for 
the pair. Thus we can understand both cases by consid- 
ering just 2 occupied bands. For example, a trivial case 
is that of two bands with no negative inversion eigenval- 
ues. This is a realization of case (i) where the inversion 
eigenvalues of each band separately multiply to +1. The 
global constraint does not prevent us from isolating all 
the occupied bands and thus, with all eigenvalues posi- 
tive, P3 = 0. All realizations of case (i) are trivial for the 
same reason. 

The first interesting case is that of two bands with a 
single pair of negative inversion eigenvalues at the same 
point (say k — 0) but positive inversion eigenvalues at 
the other point. We consider that case now. The im- 
portant consequence of the global constraint, as we just 
saw, is that the two bands with the negative eigenvalues 
can never be completely separated from each other - the 
(generically) two degeneracy points between them cannot 
be lifted or annihilated. B(k) restricted to these two- 
bands is a U (2) matrix and generically takes the form: 



where 



B(k) = e^ k \f(k)I + ig a {k)a a ) (114) 



(f(k)) 2 +g a (k)g a (k) = l. (115) 



There is a global ± sign ambiguity in the choice of f(k) 
but once the sign is chosen at one point, smoothness as- 
sures the signs at the other points. This ambiguity does 
not have implications for the final result. If we substitute 



this form into Eq. 113 all of the dependence on <p(k) (i.e. 
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the U(l) part) drops out as long as e^( fe ' is smooth (see 
Appendix [Hj Since all loops in S 3 are contractible we can 
gauge transform B(k) to remove the k-dependent phase 
so the assumption of smoothness is not an issue. What 
remains is the winding of the ST/ (2) part which must be 
an integer. Now, with only the SU(2) part we know that 
due to B(k) = St(-fc) 



f(k) = f(-k), g a (k) = -g a (-k). 



(116) 



If we think of S 3 = R 3 U {oo} then it is easy to consider 
the general form of B. The function f(k) : S 3 — > R and 
its derivative vanishes at k = (and fc = oo). This is 
a Morse function for the sphere and we can expand it 
around the origin to find: 



f(k) =N{k){M + k a C ab k b + .. 



(117) 



where C ab is a 3 x 3 constant matrix with three non-zero 
eigenvalues of the same sigrPS and N(k) is a normaliza- 
tion factor which is even in k and constrains the matrix 
B to have unit determinant. 

Without loss of generality we choose the case when 
the eigenvalues of C ab are all positive (another reason 
the eigenvalues have to have the same sign is to fix the 
boundary condition at k — oo so that it is independent 
of the path taken to get there). This choice determines 
which sign of M will lead to the non-trivial phase. Sim- 
ilarly we can expand the function 



g a (k) = N(k)(D ab k b 



(118) 



for a 3 x 3 constant matrix D ab with no restriction on 
the eigenvalues. Generically, D ab will have non-zero de- 
terminant (i.e. it will have rank 3). In cases where the 
determinant of D ab is tuned to zero, we have to look use 
a higher order Taylor-expansion in both f(k) and g a (k) - 
maintaining even terms in f(k) and odd terms in g a (k). 
As long as the boundary conditions are fixed, which re- 
quires us to keep terms in f(k) with higher order than 
g a (k), this will not change the result. Without loss of 
generality we take the case det D 7^ 0, and by rescaling 
and rotating we transform to the momentum space ba- 
sis (^1,^2,^3) with C ab = D ab = S ab . For this choice we 
have: 



N(k) 



1 



y/(M + k 2 ) 2 + k 2 



(119) 



with k 2 = k\ + k\ + k 3 . For this B(k) we have 

B(0) = sgn(A/)/ 2x2 , B(oo) = Z 2x2 . (120) 
By explicit calculation, we find that: 



(M - k 2 ) 



sign(M)-l 



Wo ((M + k 2 ) 2 + k 2 )' 

We notice that when there is an eigenvalue switch when 
passing from zero to infinity, we have P3 = 1/2, whereas 



when there is no eigenvalue switch, we have P3 = 0. Al- 
though more terms can be kept in the expansion around 
the origin this does not influence the result of the wind- 
ing number, as long as the eigenvalues of B(0) and B(oo) 
are not changed. 

We have seen from this simple argument that for two 
bands which cannot be separated, the contribution to P3 
depends only on the change in inversion eigenvalues. If 
there is only a single pair of bands which is in case (ii) 
i.e. cannot be untangled from its partner, then the other 
N — 2 occupied bands are not constrained and may each 
be isolated away from all other bands. Each of these 
isolated bands does not contribute to P3 and thus the 
non-trivial value of P3 only comes from the two tangled 
bands. To finish the proof we must consider the case 
when there are more than one set of tangled occupied 
bands. If, for example, there are four bands which have 
negative eigenvalues at k = (for simplicity we fix all the 
eigenvalues at k = 00 to be positive) we can generically 
isolate the four bands into two pairs. The two pairs can 
be separated from each other and all other bands, but 
the bands making up a single pair must share degenera- 
cies from the arguments above. Once we have decoupled 
the inversion eigenvalues of the two pairs, we can remove 
all the accidental degeneracies and isolate the pairs from 
each other because a combined pair of bands with nega- 
tive inversion eigenvalues by itself does not contradict the 
global constraint. Since the pairs can be isolated, they 
contribute independently to P3 . Each pair will contribute 
a half-integer giving P3 = n = mod Z. We see there 
is an even odd effect so that an odd number of pairs of 
bands with negative inversion eigenvalues is non-trivial 
while an even number of pairs is trivial. 

To complete the picture we will discuss how these ar- 
guments carry over to the lattice case when momentum 
space is a torus T 3 . We now have eight invariant mo- 
menta to consider which can lead to many more different 
combinations of inversion eigenvalues. We will not enu- 
merate all the possible phases but instead construct the 
necessary general principles to classify such states. We 
again consider a set of N occupied bands which does not 
have to be an even number. In general the only restric- 
tion is that the product of all the inversion eigenvalues 
of all the occupied bands is equal to one. We can gener- 
ically perform band crossings only between the occupied 
bands to split the bands into three possible classes: (i) 
a set of n+ bands with positive inversion eigenvalues at 
all ki nv , (ii) a set of n e _ bands with an even number of 
negative eigenvalues on each band, and a set of n°_ bands 
where the product of eigenvalues on each band is —1 and 
the product cannot be made equal to +1 via band cross- 
ings among the other bands in n°_. From the arguments 
above, the sets of n + and n e _ bands contribute nothing 
to P 3 , because they can be isolated one- by-one from all 
of the occupied bands. Note that the n e _ can contribute 
to non-trivial 3D QHE states in the same manner shown 
above. The number of bands in the third class n°_ must 
be an even number to satisfy the global constraint. We 
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will now show that value of P3 = (l/2)n° mod Z. We call 
this process the band decoupling process and we give ex- 
plicit examples of this band decoupling picture for lattice 
models shown in Appendix [J] 

We know that the only bands which can contribute to 
P3 are those in n°_ . We can consider these bands as a set 
of n°_/2 decoupled pairs. Each one of them contributes 
a U(2) block to the B{k) matrix. We show in Appendix 
[TJthat once the band decoupling process is finished then 
the U{\) phase in each of the n°_/2 U (2) blocks is smooth 
and can thus be completely eliminated from considera- 
tion. Thus each U(2) block can be contracted to SU(2). 
Since this implies that we are really considering maps 
from T 3 — > SU(2), which have the same dimension, we 
can connect the winding number of each SU(2) block of 
B(k) to the degree of the map. This argument follows 
along the same lines shown in Ref. 43 so we will not 
include all the details. To calculate the degree of the 
map from T 3 to an SU(2) block of B(k) we can choose 
any point in SU (2). For example, we could choose — 12x2- 
Since B(k) = B*(—k), any k which is not inversion sym- 
metric contributes to the degree of the map twice, i.e. if 
—^2x2 occurs at k it will occur at — k, thereby leaving P 3 
unchanged. The only contributions, therefore, come from 
the inversion symmetric points, and hence the winding 
number counts the number of inversion symmetric points 
that have — ^2x2 as their inversion eigenvalues. This im- 
plies that we can simply apply the Kane-Fu formula^ to 
the bands in n°_ to determine the value of P3. 



IV. SIMPLE EXAMPLE MODELS 

In this section we provide a set of explicit models that 
illustrate the majority of the technical details discussed in 
the previous sections. For each model we list the interest- 
ing phases, inversion eigenvalue structure, and describe 
what the entanglement spectra should look like. Ad- 
ditionally we provide figures showing the entanglement 
spectra for each model which confirms our analytic for- 
mulae. The models we chose are similar to ones used in 
many contexts especially in the field of topological in- 
sulators. Combined with the models introduced in Sec. 
[T] these examples provide valuable intuition about inver- 
sion invariant topological insulators and the similarities 
and differences between the states protected by inversion 
symmetry and those protected by other discrete symme- 
tries e.g. charge-conjugation, or time-reversal. 



A. ID Models 

Wc have already introduced two illustrative ID models 
for which we will analyze the entanglement spectrum. 
Additionally we will introduce a ID model of a dimerized 
chamPS 
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FIG. 5. Entanglement spectrum for the two-band ID model 
with (a) random site disorder (b) random site disorder which 
is inversion symmetric around the center of the ID chain. 
We show the entanglement spectra for many different random 
disorder configurations. 

1. Simple two band model 

Here we focus on the Hamiltonian given in Eq. [2] and 
reproduced here 

H^ik) — a cos k + sin(fc)(T 1 + (I + m — cos k)a^. 

This model has one occupied band, and in ID there are 
two inversion symmetric momenta k = 0,7r. There is a 
phase transition in this model between two insulating 
phases as a function of the parameter m and we have 
previously characterized the two phases of this model 
by examining the inversion eigenvalues of the occupied 
band. For m > the occupied band has two negative 
inversion eigenvalues. This implies that this system can 
be adiabatically continued to an atomic limit where the 
occupied atomic orbital has a negative inversion eigen- 
value. In the atomic limit the entanglement entropy is 
identically zero and so we expect on physical grounds 
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that there should be no stable entanglement eigenvalues 
at 1/2. Using our inversion criterion we see that this is 
the case since the inversion eigenvalues do not change 
between and ir. This same characterization applies for 
m < — 2 where the two occupied bands have positive in- 
version eigenvalues. The final case is for — 2 < m < 
where £(0) = —CM = +• This cannot be adiabatically 
connected to an atomic limit. Here the inversion crite- 
rion implies we should see a pair of entanglement modes 
at 1/2 which is shown in Fig. [TJ; . Thus the entanglement 
spectrum is a good indicator for a non-trivial topological 
insulator. The trivial case is shown in Fig. [TJi which has 
no 1/2 mode. We can compare these results to the values 
of the topological invariant x which takes on values x — 
when to < —2 or to > and x = 1 when —2 < to < 0. 
The number of required 1/2 modes in the entanglement 
spectrum is 0, and 2 respectively which is exactly what 
we found numerically. 

To illustrate the protection due to inversion symmetry 
we also consider H 2 with an onsite disorder term added 



H = Ho 



(122) 



where the Wi are randomly chosen from a uniform distri- 
bution [— W/2, W/2}. If it is purely random, uncorrelated 
disorder the mid-gap entanglement modes are no longer 
protected as shown in Fig. [5^,. Next we mirror the dis- 
order around the center of the chain to make a system 
with inversion symmetric disorder. Although unphysical, 
this helps illustrate the fact that only inversion symme- 
try is required for the protected mid-gap entanglement 
states as we show in Fig. [Hfj where the cut is along the 
remaining inversion center. 



(b) 




FIG. 6. (a) Id Dimerized Chain (b)2d Dimerized square lat- 
tice model. Solid and dotted lines represent different hopping 
amplitudes. 
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FIG. 7. Entanglement Spectra for (a)Dimerized 

Chain(b)Dimerized Square Lattice 



2. Simple four band model 

The Hamiltonian for the simple four band model was 
given in Eq. [5] and the set of inversion eigenvalues for the 
different phases were listed in Eqs. 6]|l0 For convenience 
we reproduce the Hamiltonian here 

H 4 (k) = sin(fc)f 1 + sin(fc)f 2 
+ (1 — to — cos/s)f + ST 2 4 + ecos(fc)(l + f ). 

From the inversion criterion we see that Case 2) and Case 
3) should have entanglement modes at 1/2. In fact, Case 
2) should have a pair of modes localized on each cut, one 
for each occupied band that flips the sign of the inversion 
eigenvalue. The entanglement spectra for the five cases 
are shown in Fig. [2j-h and they agree with the analytic 
prediction. For cases 1-5 the invariant x = 0, 2, 1, 0, 
yielding 0,4,2,0,0 entanglement modes at 1/2 which is 
what we found numerically. 



3. Dimerized Chain 

As a final ID test case we will look at spinless fermions 
hopping on a dimerized chain. This model is the familiar 
Su-Schrieffer-Heeger model for electrons in a polyacety- 
lene chainP^ Later we will extend this model into 2D and 
3D to illustrate anisotropic systems with non-trivial en- 
tanglement spectra. The layout of the chain is shown in 
Fig. [6^i along with the choice of two-atom unit cell. The 
Hamiltonian is given by 



H = J2 A \ c La c ™a - c\ nB c mB 

rn 

+ ((-* - <5)<4 A c mB + (-t + b)c\ iB c m+ 



h.c. 



(123) 



where A, B indicate sublattice A or B and t > 0. For 
our purposes we will set the onsite energy A = 0. With 
this choice the Hamiltonian has an inversion symmetry 
with respect to the middle of a bond with V — <j x i.e. V 
exchanges sublattices A and B. 
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The Hamiltonian can be Fourier transformed and be- 
comes 

tf=5>U(-(* +*)-(* 



+ (t - S) sin ka v ] 

*fe = (c k A C kB ) T ■ 



(124) 

This model has two insulating phases {1)5 < (ii) S > 0. 
At the two inversion invariant points we have the Bloch 
Hamiltonians H(0) = —2ta x and H(tt) = —2Sa x . As ex- 
pected they both commute with V . For fixed t = 1 we 
see that the inversion eigenvalue of the occupied band 
at fc = is fixed to be +. For fc = it the eigenvalue 
depends on the sign of 5. So for 5 < 0(> 0) we have 
C(7r) = — 1(+1). We see that the non-trivial phase occurs 
when 5 < 0. In this case the Wannier center of the elec- 
trons is shifted to the mid-bond center between unit cells. 
Thus if we cut the system between unit cells there will 
be a charge polarization. If 6 > the Wannier center is 
on the mid-bond center within a unit cell. Thus our def- 
inition of unit cells is important and simply reflects the 
fact that the polarization is not well-defined absolutely 
but is gauge- variant. 

The entanglement spectra for this model for 5 < is 
shown in Fig. [7^,. We expect that because the inversion 
eigenvalues change from positive to negative between fc = 
0, and 7r that there should be 1/2 modes in the spectrum 
and this is confirmed numerically as seen in the figure. 
The topological invariant x = 0, 1 for <5 > and 5 < 
respectively. The number of the entanglement modes in 
the spectrum matches the value of 2% as expected for a 
system with periodic boundary conditions and thus two 
entanglement cuts. 

There is a subtlety in the entanglement characteriza- 
tion of this model which we will now discuss. The real- 
space Hamiltonian as written in Eq. 123 has broken 
translational symmetry. As stated earlier, our classifi- 
cation method only applies to translationally invariant 
Hamiltonians since we need to evaluate inversion eigen- 
values in momentum-space. However, by construction, 
the Hamiltonian for the dimerized chain only has a very 
mild breaking of translational symmetry. In fact, as im- 
plicitly assumed in the analysis, we can just define a unit 
cell encapsulating two-sites, and in terms of the larger 
unit cell the Hamiltonian is translationally invariant and 
our method applies. Our choice for a unit cell has al- 
ready been implicitly assumed by the time we write the 
Bloch Hamiltonian in Eq. |124| The choice made here is 
that sites connected by the hopping term — (t + S) form 
a single unit cell and the hopping between unit cells is 
(—t + S). After making this choice we are free to carry 
out the entanglement analysis by cutting the system be- 
tween unit cells. Cutting the system within a unit cell 
is not a position-space cut of a translationally invariant 
Hamiltonian. 

Now, the important outstanding question is: is the 
analysis invariant under the choice of unit cell? It is, and 
we explicitly show it for this model. Suppose that we 



take the unit cell to be sites connected by the hopping 
(—t + S). Then the Bloch Hamiltonian becomes 

H = E*! [H*-*) -(* + <*) cos k)a x 



+ (t + S) sinfc a y ] <3> fe 
H(0) = ~2ta x H(tt) = 2Sa x . 



(125) 



If we fix t > then the inversion eigenvalue for the 
occupied band at fc = is always positive. Thus we 
see that for 5 > the system is in a non-trivial phase 
with opposite inversion eigenvalues at and 7r, and for 
S < the system is trivial. The sign of S that exhibits 
the non-trivial phase has changed when compared with 
the choice of the other unit cell. The physics, however, 
remains identical. Now when S > the Wannier cen- 
ters are shifted to the bonds between the new unit cells. 
This would have lied within the unit cell for our previous 
choice and explains why the sign of S has changed. For 
the choice of unit cell in H the entanglement spectrum 
will have mid-gap modes when S > i.e. when the in- 
version eigenvalues are opposite at the two invariant mo- 
menta. Thus we see that for the new choice of unit cell 
the physics and entanglement analysis yields the same 
results. 



B. 2D Models 



1. Dimerized Square Lattice Model 

The first 2D model we consider is a trivial extension 
of the dimerized chain as illustrated in Fig. [6}d. This 
extension has a Hamiltonian which is simply constructed 
from Eq. [124] 

H = E *I [(-0 + S)-(t-S) cos k x ) a x 

k 

+ (t - S)smk x a y - 2t y cosky] * fc . (126) 

This model has an inversion symmetry with V = a x . At 
the inversion invariant points we have 

H(0,0) = -2t y - 2ta x 
H(tt,0) = -2t y - 2Sa x 
H(0,tt) = 2t y - 2ta x 
#0,7!-) = 2t y - 2Sa x . 



We see that although this is more complicated than the 
ID case everything still commutes with V . For simplic- 
ity we pick t — 2t y — 1 and focus on the two gapped 
phases 8 < and S > 0. In these two phases we have the 
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following set of eigenvalues 



Hamiltonian is 




H = 



6>Q: 



S < : 



The product of all the eigenvalues in each case is +1 
so the parity of the first Chern number for these two 
cases is even. In fact, it is exactly zero for this model. 
Just as in the ID case we see that the S < phase is 
interesting. Here the Wannier center for each electron 
is moved along the x-axis to the mid-bond center be- 
tween each unit cell. If we cut an edge which is per- 
pendicular to the x-direction there will be a finite charge 
density on the boundary. Although the inversion sym- 
metry is not enough to determine the polarization we 
see that this Hamiltonian also has a reflection symme- 
try MH(k x , fc^VW- 1 = H(-k x , k y ) with M = a x = V. 
Thus, the charge polarization in the x-direction is quan- 
tized and equal to P\ — e/2a where a is the lattice con- 
stant in the y-direction. 

Finally we can consider the entanglement spectra of 
these two phases. For 5 > the phase is adiabatically 
connected to the atomic limit and thus will not require 
the existence of 1/2 modes. For 6 < we must first 
specify a cut direction to locate the 1/2 modes. Let us 
first pick the cut to be parallel to the x-direction. Thus k x 
remains a good quantum number and we can ask whether 
or not there are 1/2 modes at k x — or k x = tt. For 
k x = we look at £(0, 0)C(0, tt) = +1 and for k x = tt 
we consider ((n, 0)((ir 7 tt) = +1. Thus for this cut there 
are no 1/2 modes. Next we look at a cut parallel to the 
y-direction such that k y is a conserved quantum number. 
For k y = we have C(0, 0)£(7r, 0) = —1 and for k y = tt we 
have £(0, tt)£(tt, tt) = —1 which implies there will be 1/2 
modes at both k y = and tt. The entanglement spectrum 
for a cut parallel to the y-direction is shown in Fig. [7]d. 
In this figure there are clear 1/2 modes at k y = 0,tt. 
In fact, for this simple model there are 1/2 modes for 
all values of k y though our criterion only constrains the 
modes at the inversion invariant points. 



2. Chern Insulator 

Next we move on to study the well-known 2D topologi- 
cal insulators beginning with the Chern insulator (quan- 
tum anomalous Hall effect) P This is a 2D topological 
insulator which exhibits a quantum Hall effect and is 
classified by an integer invariant, the Chern number.^ 
Instead of studying the initially proposed honeycomb lat- 
tice model we will use the square lattice version which is 
nothing but massive Dirac fermions on a lattice. The 
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m+l,n 



(ia x - a z )* T 



(2-m)^^a 2 * m ,„} 
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We can Fourier transform to get the Bloch Hamiltonian 
H(k) = sin k x a x + sin k y a y + M(k)a z (128) 



M(k) = 2 — to — cos k x — cos k y 



(129) 



This model exhibits several different phases as a function 
of m. For to < 0, to > 4 the system is in a trivial insu- 
lator phase and for0<m<2,2<m<4 the system 
is in Chern insulator (quantum Hall) phases with Chern 
number —1 and +1 respectively. This Hamiltonian has 
an inversion symmetry with V = a z and at the four in- 
version invariant momenta we have 

#(0,0) = -mo z 
H(n,0) = 2-ma z 
H(0,tt) = 2-ma z 
H(n, 7r) = 4 - ma z . 

For to < 0, to > 4 the inversion eigenvalues are all posi- 
tive/negative respectively. For < to < 2 C(0, 0) = — 1 
and C( 7r I 0) = C(0, tt) — ({tt,tt) = +1 and the system 
must have a Chern number with odd parity. It does 
since C\ = —1. For 2 < m < 4 all the eigenvalues except 
C(7r,7r) are negative and again the Chern number must 
have odd parity (C\ = +1). 

The location of the 1/2 modes in the entanglement 
spectrum is also clear. If we choose to cut along the x 
or y directions the picture remains the same. This indi- 
cates that we are not dealing with an (weak) anisotropic 
insulator as in the dimerized case but a fully 2D topolog- 
ical insulator state. This is similar to saying that for the 
quantum Hall effect, no matter what edge we cut in the 
system, there will be edge states. For definiteness assume 
that we cut parallel to y so that k y is a good quantum 
number. For < 777 < 2 there will be a 1/2 entanglement 
mode at k y = and for 2 < to < 4 there will be a 1/2 
mode at k y = tt. In addition to these modes there is ac- 
tually a dispersing set of modes that are localized on the 
cut. To clearly see the dispersing modes we look at the 
entanglement 'energies' which are defined to be 
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where £ m are the eigenvalues of the reduced correla- 
tion matrix Cl. The entanglement energies show the 
full structure of the entanglement spectrum because they 
clearly separate the eigenvalues of Cl which are clustered 
near zero and one. The energy, entanglement eigenval- 
ues, and entanglement energies for the Chern insulator in 
the to<0, 0<to<2 and 2 < to < 4 phases are shown 
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FIG. 8. Energy spectrum with an open boundary, entanglement spectrum with a cut parallel to the y-direction, and entangle- 
ment energies for (Left panel a,d,g)A trivial insulator (Middle panel b,e,h)Non-trivial Chern insulator with negative inversion 
eigenvalue at (k x , k y ) — (0, 0) implying entanglement mid-gap modes at k y — (Right Panel c,f,i)Chern insulator with negative 
inversion eigenvalues at (k x ,k y ) — (0,0), (it, 0) and (0, 7r). This implies entanglement mid-gap modes at k y — ir. 



in Fig. [8] In the two non-trivial phases the location of 
the 1/2 mode is different. For < m < 2 it is at k y = 
and for 2 < m < 4 it is at k y = it. For the entanglement 
energies these become zero modes. These entanglement 
spectra were cut from a torus geometry so that there are 
two entanglement cuts. This is the reason why there are 
entanglement modes dispersing in both directions in the 
C\ 7^ phases. 



3. Quantum Spin Hall Insulator 

The quantum spin Hall insulator (QSH ) is a time- 
reversal invariant topological insulato j 2 * 4 * 5 * which is most 
easily thought of as two copies of the Chern insulator, 
one for each spin, with opposite chiralities. The realis- 
tic material in which this state is realized, HgTe/CdTe 
quantum wells, is best modeled by exactly this type of 
Hamiltonian. The effective HgTe Hamiltonian is a four- 
band model on the square lattice with a Hamiltonian 



given by 

H(k) = sin fea.fi + sin k y T 2 + M (k)t (131) 
M(k) = 2 — m — cos k x — cos k y . (132) 

where f i = a z ®t x ,T 2 — 1 <E> t v , f = 1 ® t z where a a is 
spin and r a is the orbital degree of freedom. For this sys- 
tem the time-reversal operator is T = (ia v ' ®1)K and the 
inversion operator is V = f q. This Hamiltonian is invari- 
ant under both symmetries. It exhibits phases directly 
analogous to the Chern insulator i.e. it is a trivial insu- 
lator for m < 0, m > 4 and a topological quantum spin 
Hall insulator for < m < 2 and 2 < m < 4. The only 
difference with the Chern insulator is that now there are 
two occupied bands which are related by time-reversal 
symmetry. The total Chern number of the ground state 
vanishes but there is a Z 2 invariant given (in the presence 
of time-reversal and inversion^!) by 

Xz 2 - 1] Cn(0,0)C„(7r,0)C n (0 J 7r)C„(7r J 7r)(133) 

n£occ./2 

where the product runs over half the occupied bands. 

(2) 

This invariant has the same formula as xp defined 




FIG. 9. Energy spectrum with an open boundary, entanglement spectrum with a cut parallel to the y-direction, and en- 
tanglement energies for (Left panel a,d,g)A trivial insulator (Middle panel b,e,h)Non-trivial quantum spin Hall insulator with 
time-reversal symmetry (Right Panel c,f,i)Quantum spin Hall insulator with mildly broken time-reversal symmetry. Comparing 
h and i one can see that all Kramers' degeneracies are lifted when time-reversal is broken except the ones at e = 0. By just 
looking at e and f it is difficult to tell the difference in the two cases i.e. that spectral flow is broken. 



above, but we distinguish it here to prevent confusion 
since this invariant only implies a non-trivial physical re- 
sponse when time-reversal is preserved. To specify which 
half of the occupied bands you just take one from ev- 
ery Kramers' pair of bands. If we focus on the transi- 
tion when m ~ we see that for m < the inversion 
eigenvalues of both occupied bands are all positive. For 
2 > m > the inversion eigenvalues at (k x ,k y ) = (0, 0) 
for both occupied bands are negative while all others are 
positive. The product over inversion eigenvalues of all 
the occupied bands is trivial, but if we only multiply 
over half the Kramers' pairs we find that \z 2 — ~ 1 an d 
is non-trivial. Since the product over all the bands is triv- 
ial this means that the parity of the first Chern number is 
even, in fact it is zero for this case. The energies, entan- 
glement eigenvalues, and entanglement energies for these 
two phases are shown in the left and middle panels of 
Fig. [9j Without time-reversal symmetry this type of in- 
version invariant (product over half the occupied states) 
does not uniquely specify a topological response in 2D 
since a C\ = 2 quantum anomalous Hall state could have 



the same inversion eigenvalue structure. 



4- Quantum Spin Hall Insulator without time-reversal 
symmetry 

So far the studies of the Chern Insulator and QSH 
insulator have just been reconfirmed by recognizing the 
importance of inversion eigenvalues when there is an in- 
version symmetry. The most interesting prospect is when 
we take the QSH effect and break time-reversal but keep 
inversion. The importance of inversion symmetry for this 
type of case in 3D was emphasized in Ref. [29] To break 
time-reversal symmetry we will consider an additional 
Zeeman term in the QSH Hamiltoniarn^ 




sin k y T 2 



M(k)T 



(134) 
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For small values of B x this term lifts the Kramers' de- 
generacy of the occupied bands but does not cause any 
crossings at the Fermi-level. Although time-reversal is 
broken, inversion is still preserved and we can still see 
that xz 2 = ~ 1- This is still well defined because the 
product of all inversion eigenvalues at a particular k% nv is 
still trivial for all k inv . Thus, this system is an inversion 
invariant topological insulator. It was first noted that 
such states could exist in Ref. [2^1 where it was suggested 
that as long as inversion symmetry was not broken the 
entanglement spectra for the time-reversal preserved and 
broken cases were the same. However, there is an impor- 
tant difference between the two. For the time-reversal 
broken QSH state we show the entanglement eigenvalues 
and entanglement energies in the right panel of Fig. [9] 
Comparing with the middle panel we see that the entan- 
glement eigenvalues seemingly show very little difference 
due to the fact that most are exponentially close to zero 
or one, but the entanglement energies are quite different. 
The mode at 1/2 is protected by inversion symmetry but 
all of the other "Kramers'" degeneracies are lifted e.g. all 



the crossings at k y 



-7r and k,, = ir are lifted. This oc- 



curs because the spectral flow between the bulk valence 
and conduction bands is cut-off when time-reversal is bro- 
ken. The edge states no longer tie together both bands 
and there is no "anomaly" -type structure. Thus states 
on the left and right half of the system are no longer tied 
together through the bulk in a topological way. 

We ask now if there is anything interesting in this sys- 
tem once time-reversal is broken? As long as inversion is 
preserved we cannot connect this state to a trivial atomic 
limit while preserving inversion symmetry and there must 
always be a finite entanglement entropy. The finite en- 
tanglement entropy is due to the fact that the mode at 
1 /2 is protected and cannot be removed. Thus the system 
is not trivial in the sense that it cannot ever be continu- 
ously deformed to a trivial atomic insulator; the entangle- 
ment spectrum clearly shows this. Now we can also ask 
if there is any non-trivial physical response? The robust 
electromagnetic response discussed in Refs. HH1 and 
[50] comes from coupling the quantum spin Hall state to 
varying adiabatic parameters. For example, applying a 
magnetic domain wall to the edge of a QSH system in- 
duces a mid-gap state and fractional charge localized on 
the domain wall. This response requires time-reversal 
symmetry to be broken and thus could still exist in the 
presence of a Zeeman term. The Zeeman term is an ex- 
ternal time-reversal breaking field and it opens a gap in 
the edge states. If we are able to create a magnetic do- 
main wall on the edge, i.e. a field strong enough to re- 
verse the direction of the Zeeman field in some region of 
the edge then there will be trapped domain wall states. 
Thus the state is a topological insulator in a very phys- 
ical sense as well. Unfortunately the inversion invariant 
{x-p ) m 2D does not imply that this must be the phys- 
ical response. As mentioned above a C\ = 2 quantum 
anomalous Hall effect can have the same value of the Z 2 
invariant. However, for this model Hamiltonian we are 




FIG. 10. Energy spectrum with open boundary conditions 
for (a)8-band model with T symmetry (b) 8-band model 
with random inversion preserving perturbation. Entangle- 
ment spectrum for (c) 8-band with T (d) 8-band with random 
inversion preserving perturbation. 



saved because there is an additional mirror symmetry 
MH(k x ,k y )M^ 1 = H(k x ,-k y ) with M = a x ®t z . Note 
that [M,V] = [M,T B ] = 0. Thus, the inversion eigen- 
values are still valid labels and the Zeeman term does 
not break the mirror symmetry. This symmetry forbids 
a non-zero C\ and thus the QSH response is the unique 
result. 



5. 2D 8-band model 

Here we consider a more complicated case of a model 
with 8-bands total and four occupied bands. The model 
we choose is simply two-copies of the QSH model. The 
Bloch Hamiltonian is given by 

H(k) = sin k x j x + sin k v ~f v + M(k)Y (135) 



with M(k) given in Eq. 132 and 7 



This model preserves time-reversal symmetry with T = 
(1 <g> ia y <S> 1)K and inversion symmetry with V — Y~ ■ in 
the presence of time-reversal symmetry this model does 
not yield any non-trivial topological insulators since you 
always get an even number of pairs of edge states. If one 
adds perturbations which break time-reversal symmetry 
then it is possible to generate non-trivial states such as 
Chern insulator states. At half-filling this model will 
have four occupied bands, and will exhibit edge states 
for the same range of parameters as the quantum spin 
Hall model above in Eq. 131 These edge states are not 



protected generically i.e. one can add a time-reversal and 
inversion preserving perturbation to the model which will 
open a gap in the edge states. We use this model to il- 
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FIG. 11. (a)Energy spectrum for the small anisotropic 8-band 
Hamiltonian with periodic boundary conditions. Difference in 
negative parity eigenvalues at k y = 0,-k is three, (b) Entan- 
glement spectrum for the same. At k y = there are six modes 
at £ = 1/2. The lines are only filled in as guides to the eye. 
The only allowed k- values are k y — and k y — it. This is why 
the spectra appear to have kinks. 



lustrate three points (i) even though there is no topolog- 
ical invariant in the system associated with time-reversal 
symmetry, the presence of inversion symmetry predicts 
that there will be non-trivial mid-gap states in the en- 
tanglement spectrum (ii) these mid-gap states are com- 
pletely stable as long as you do not break inversion sym- 
metry (iii) the number of mid-gap states is proportional 
to the difference in negative parity eigenvalues at the in- 
variant momenta. 

In Fig. [lO] we show the energy and entanglement spec- 
tra for two different cases. In Fig. [I0"fc,,c we simply diag- 
onalizc Eq. |135| with m = 1.0 on an open boundary. One 
can clearly see the (unstable) edge states lying in the 
mid-gap region. The entanglement spectrum is shown 
for a system with periodic boundary conditions cut par- 
allel to the y-axis. The difference between the number 



of negative inversion eigenvalues at (k x ,k y ) = (0,0) and 
(tt, 0) is four giving a total number of 2 x 4 = 8 modes 
at 1/2. If we leave out the identity matrix, and V it- 
self, there are 30 matrices which commute with V . To 
illustrate the stability of the 1/2 modes in the entangle- 
ment spectrum we add a perturbation which includes all 
30 matrices with random couplings chosen from a uni- 
form distribution [—6/2, 8/2] where 6 > is chosen small 
enough not to close the bulk gap. This perturbation will 
break all of the 'accidental' symmetries in the problem 
but preserves inversion symmetry. We see in Fig. |10| 
where 6 = 0.20 that the energy spectrum still has states 
in the gap but the crossing-points have been lifted. The 
entanglement spectrum, however, still has eight exact 1/2 
modes. 

Finally we take a very small anisotropic sized system 
with L x — 20 and L y — 2. We cut the system in the 
middle between x = 10 and x = 11 and plot the en- 
tanglement spectrum vs k y . There are only two allowed 
values for k y = 0, n. The energy spectrum for such a sys- 
tem (with the random perturbation included but chosen 
from a uniform distribution [0, 6] which is no longer sym- 
metric about zero) is gapped and has an entanglement 
spectrum with six mid-gap 1/2 modes when 8 = 0.19. 
Counting the occupied states we find numerically that 
there is a difference of three negative inversion eigenval- 
ues which agrees with the entanglement spectrum. The 
energy and entanglement spectra are shown in Fig. 11 



C. 3D Models 

1. Dimerized Cubic Lattice 

The 3D dimerized model on a cubic lattice is a trivial 
extension of the 2D dimerized case into 3D. The Hamil- 
tonian is given by 

H = ^2^\ [(-(* + <*) -(t-6) cos k x )a x 

k 

+ (t — 6) sin k x a y — 2t y cos k y — 2t z cos k z ] ■ 

(136) 

For t — 2t y — 2t z — 1 there are two different phases 
8 < and 6 > 0. As before, for 6 > the Wan- 
nier center of the electron is located within a unit cell 
and all inversion eigenvalues are +1. For 8 < the 
Wannier center is shifted along the x-axis to the mid- 
bond site between unit cells. The inversion eigenval- 
ues are C(000) = C(OOtt) = C(OttO) = ((Otttt) = +1 
and C(ttOO) = C^Ott) = C^ttO) = ((tttttt) = -1. 
Again this system has a non-trivial charge polarization 
on a surface with X CIS 3b normal vector. However, this 
is protected by the reflection symmetry about the yz- 
plane i.e. MH(k x , k y , k z )M~ x = H(—k x .k y .k z ) with 
M = cr x = V . Thus, since the product of the reflec- 
tion eigenvalues is — 1 for each reflection invariant line in 
the Brillouin zone, the polarization on a surface perpen- 
dicular to the x-axis is P\ = e/2a 2 where a 2 is the area 
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FIG. 12. Entanglement Spectra for (a)3DQHE (b)3d WTI 
with a cut parallel to the x-y plane plotted along a line in the 
Brillouin zone. The states in (b) are all doubly degenerate 
compared to figure (a). 



of a plaquette in the yz plane. The product of all the 
parity eigenvalues is trivial as it must be in 3D, and ad- 
ditionally the product of the eigenvalues in every plane is 
trivial. For an entanglement cut such that k y and k z are 
good quantum numbers it is clear that there will be 1 /2 
modes at (k y ,k z ) — (0, 0), (tt, 0), (0, 7r), and (7r, tt). On 
the other translationally invariant cuts parallel to the xz 
or xy planes there will be no 1/2 modes. 



2. 3D Quantum Hall Effect 

The 3D quantum Hall effect state can be thought of 
as stacks of 2d quantum Hall states which are connected 
together. We will use a very simple model for the 3D 
quantum (anomalous) Hall effect which is a trivial ex- 
tension of the 2D Chern insulator. The Hamiltonian is 



H(k) = smk x a x + smk y a y + M(k)a z 



(137) 



M(k) = 2 — in — cos k x — cos k y — t± cos k z . (138) 

This system has an inversion symmetry with V = u z . 
This model exhibits several different phase transitions 
but we will only focus on one, namely the phase transition 
that occurs with a gapless point at [k x , k v , k z ) = (0, 0, 0). 
For m < — tj_ the system is in a trivial insulating phase 
with all inversion eigenvalues positive. At m = — 1± the 
system becomes gapless and stays gapless until m > tj_- 
For m > t± the system is in a 3D quantum Hall effect 
phase^D with effectively 2D quantum Hall states stacked 
up in the z-direction. The Bloch Hamiltonians at the 
inversion invariant points which switch eigenvalues are 

#(0,0,0) = -(to-M_lK 
#(0,0, tt) = {-m + ti_)a z 

At to = — tj_ the eigenvalue around the T-point switches 
from positive to negative but the system is still gapless. 
Then at m = +t± the eigenvalue at (0, 0, tt) switches and 
the system becomes a gapped insulator. The product 
over all the eigenvalues is trivial as expected but if we 



restrict the product to the k z = or k z = it planes the 
product is negative. We proved earlier that this indicates 
a non-trivial 3D QHE response. 

For this arrangement of eigenvalues we can calculate 
the location of the 1/2 modes in the entanglement spec- 
trum. For a cut parallel to the a;j/-plane there will be no 
1/2 modes at the inversion invariant points. If we take 
the cut e.g. parallel to the zx-plane then there will be 
1/2 modes at (k z ,k x ) = (0,0) and (tt, 0). These nodes 
are shown in Fig. [12^ ,. The figure shows, in addition to 
the 1/2 modes at the two invariant momenta, a line of 
1/2 modes between (k z ,k z ) = (0,0) and (tt, 0). 



3. 3D Weak Topological Insulator- 
There are several different classes of the recently pro- 
posed 3D time-reversal invariant topological insulators. 
The anisotropic classes, the so-called weak topological 
insulators, are effectively 2D quantum spin Hall states 
stacked into 3D. This is similar to the 3D quantum Hall 
effect and is essentially just two copies of that system, 
one for each spin. We use the following model: 

H(k) = sin k x t ! + sin k y t 2 + M (fc)f (139) 
M[k) = 2 — 777 — cos k x — cos k y — t± cos k z . (140) 

where the T a matrices are the same as in the quantum 
spin Hall state. This system is time-reversal and inver- 
sion invariant with an inversion operator V = Tq. It 
exhibits phase transitions at the same values of m as 
the 3D quantum Hall effect and the only difference is 
that there are two occupied bands instead of one. For 
2 — tj_ > 777 > t± The system has two pairs of nega- 
tive inversion eigenvalues, one pair at (0, 0, 0) and one 
at (0,0, tt). The rest of the eigenvalues are all positive. 
The total product of inversion eigenvalues is trivial, and 
unlike the 3D quantum Hall case the product of the 
eigenvalues when restricted to the k z — 0, tt planes is 
also trivial. However, there is still something non-trivial 
here which arises from taking the eigenvalues from only 
one of the Kramers' pairs at each invariant momentum. 
We see that this product is non-trivial and indicates an 
anisotropic inversion invariant topological insulator. In 
this case, since time-reversal symmetry is preserved, it 
is a weak topological insulator stateP However, with- 
out time-reversal symmetry the inversion invariant be- 
ing non-trivial does not require it is a weak topological 
insulator. As a counter example it could be a 3D quan- 
tum Hall effect with a Chern number "per layer" that 
is an odd multiple of two (unless there is a reflection 
symmetry which requires the Chern number to vanish 
in each plane). The interesting thing about this model is 
that even when time-reversal is softly broken (i.e. broken 
without causing a phase transition) the system still is not 
in a trivial topological state and even though the surface 
states are no longer protected it can exhibit non-trivial 
behavior in the entanglement. The entanglement spec- 
trum for the T-invariant case is shown in Fig. |12b and 
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FIG. 13. (Left column) Strong topological insulator (a) 
Energy spectrum plotted along a line in the Brillouin 
zone for open boundary conditions along the z-direction 
(c)entanglement spectrum with a cut parallel to the x-y plane 
with periodic boundary conditions (e)entanglement energies. 
(Right column) Strong topological insulator with time rever- 
sal breaking (b)energy spectrum (d)entanglement spectrum 
(f) entanglement energies. In (a) and (b) the mid-gap states 
are localized on the surfaces. In (b) there is a gap in the sur- 
face states due to time-reversal symmetry breaking. In (f) all 
Kramers' degeneracies are lifted except the one at e = 0. 



exhibits the same 1/2 mode structure as the 3D quantum 
Hall effect model but with twice as many modes. 



4- 3D Strong Topological Insulator 

The last class of models we will consider is the 3D 
lattice Dirac model which is the minimal model for time- 
reversal invariant strong topological insulators in 3D. The 
Bloch Hamiltonian is given by 



H (k) = sin k x T\ 
M (k) = 3 - m - 



+ sin k y f 2 
cos k x — cos k y 



smfe 2 f 3 + M(fc)f( 
cos k z 



(141) 



where f 3 = a v (g) t x . As a function of m this model ex- 
hibits many phase transitions. We will focus on the range 
m < and < m < 3. There is a phase transition at 
m = with a band crossing at the T-point in k-space. 
Four bands meet at this point and a pair of inversion 
eigenvalues is exchanged. For m < the inversion eigen- 
values are positive for both occupied bands at all the 
invariant momenta. The Bloch Hamiltonian at k = is 
H (0) = — mTo and thus when m switches sign the inver- 
sion eigenvalues at k — are exchanged. For < m < 3 



both inversion eigenvalues at k = are negative. In 
this phase the product over all inversion eigenvalues is 
trivial, but if we only keep one of the Kramers' pairs, 
then the product of all the eigenvalues of half the occu- 
pied bands is non-trivial. In the presence of inversion 
and time-reversal, which is the case here, this invariant 
is the strong topological Z 2 index. 7 Physically this index 
has two implications: (1) the presence of an odd- number 
of massless Dirac cones on any surface (2) a topological 
magneto-electric effect. To see the topological response 
one must apply a time-reversal breaking field on the sur- 
face to open a gap in the gapless Dirac fermions. This 
induces a quantum Hall effect confined to the surface 
which leads to the magneto-electric response. For the 
topological phase we picked (0 < m < 3), the surface 
states are located around the surface T-point and the 
(pair of) entanglement modes will be located at the T- 
point of the conserved momenta parallel to the entangle- 
ment cut. The energy spectrum, entanglement eigenval- 



ues, and entanglement energies are shown in Fig. 13 j,,c,e 
respectively. The energy spectrum is shown with periodic 
boundary conditions along x and y and open boundary 
conditions along z. The surface states are shown in red 
and the spectrum is plotted along a ID path in the Bril- 
louin zone. 



5. 3D Strong Topological Insulator without time-reversal 
symmetry 

Now we can consider the more interesting question of 
the properties of the system when we break time-reversal 
but keep inversion. We add the same Zeeman term 
shown in Eq. |135| to the bulk of the insulator. We only 
break time-reversal softly which opens a gap in the sur- 
face states but does not close the bulk gap. Thus, while 
we can no longer consider the Z 2 invariant protected by 
time-reversal, this system will still exhibit a magneto- 
electric effect since the Zeeman field simply establishes 
a quantum Hall effect on the surface. This distinction 
was first considered in Ref. Additionally, although 
we can add surface potentials to push the surface states 
into the bulk bands we still cannot adiabatically con- 
nect this insulator with an atomic limit. This is clearly 
shown in the entanglement spectrum where there are still 
modes protected at 1/2 which cannot be removed with- 
out breaking inversion or passing through a phase transi- 
tion. Thus, it is impossible to continuously connect this 
to a system with vanishing entanglement entropy. Fig. 
13p,d,f show the energy, entanglement eigenvalues, and 
entanglement energies respectively. For the energies the 
red surface states are now clearly gapped, but the en- 
tanglement eigenvalues are hard to distinguish between 
the time-reversal invariant and breaking cases. The dis- 
tinction, we see, comes when we look at the entangle- 
ment energies which show that while there is still a pair 
of zero modes, all of the other degeneracies at the in- 
version invariant momenta (which arose from Kramers' 
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degeneracies) are lifted generically. This shows that the 
spectral flow between the valence and conduction bands 
of the entanglement spectrum has been cut off. 

Now we will prove that P3 7^ for this model. To be 
explicit we take 

H (k) = sin feiTi + sin ^2^2 + sin k^T^ 

. (142) 
+(5/2 — cos k\ — cos faz — cos A^To + ^Tb 

This Hamiltonian can be connected through a gapped 
interpolation with the trivial Hamiltonian: 

H (k) = sinfcif! + sin fc 2 f 2 + sin fc 3 f 3 + 5/2f (143) 

using the inversion symmetric homotopy 

H(k,0) = |(1 + cos9)H(k) 

2 (144) 
+ |(1 - cos 6) H (k) + sin6»f 5 

where I5 = r r 1 r 2 r 3 . The second Chern number gen- 
erated by H{k,&) is C2=l (odd) and consequently P 3 — 
I mod(Z) for (A) since we specifically chose the form of 
Ho(k) to be trivial. Note, that as with the cases shown in 
Appendix [j] we calculated C2 numerically using the stan- 
dard gauge invariant formula in terms of ground-state 
projection operators. 

V. CONCLUSIONS 

The question of what makes an insulator "topologi- 
cal" has many answers. In this article we presented an 
answer which encompasses all of the known topological 
insulators. The fundamental distinction between an or- 
dinary band insulator and a topological insulator is the 
inability to adiabatically connect a topological insulator 
to the atomic limit. This distinction can have many man- 
ifestations including non-trivial topological responses to 
external fields and robust boundary states, however these 
properties are not necessary conditions for a topological 
insulator. In fact, we have seen examples in this paper 
without protected boundary states, and examples with 
no topological response. At first sight these insulators 
seem to have no characteristics which distinguish them 
from trivial band insulators. Admittedly these are not 
the most interesting systems to consider experimentally, 
but they still show a striking signature in the entangle- 
ment spectrum. In fact, all known topological insulators 



show a signature in the entanglement spectrum when the 
bi-partition is a position-space cut. Instead of taking the 
whole spectrum, one can calculate just the entanglement 
entropy, which for all topological insulators cannot be 
adiabatically deformed to zero. This fact is what serves 
as the basis for our definition and unifies the inversion 
symmetric insulators with the ones that are invariant 
under other discrete symmetries. The experimental rele- 
vance of inversion symmetric topological insulators is un- 
clear, but not out of the question.^ Although in principle 
disorder immediately destroys any stability of the topo- 
logical state (unlike the typical topological insulatorsES) , 
the robustness of the insulator state is ultimately a ques- 
tion to be answered in practice. Many of the well-known 
topological insulators simplify when inversion symmetry 
is required along with the discrete symmetry that stabi- 
lizes the topological state. Thus it seems like the most 
interesting inversion symmetric insulators are ones which 
are derived from parent topological insulator states with 
weakly broken T or C symmetries. These types of mate- 
rials would be the first place to search for signatures of 
topological protection due to inversion symmetry. 

Note: During the last stages of preparation of this 
manuscript, we became aware of a paper by A. Turner, Yi 
Zhang, R. Mong and A. Vishwanath dealing with similar 
issues. 
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Appendix A: Proof that the one-body correlation 
function is a projection operator 

Here we prove that the one-body correlation function 
over the full system (not only over part of the system) is 
a projector. This can easily be proved: Cy (which is a 
matrix at each i,j i.e. C 1 "- ) has the property: 



Ei C l3 {ky)C jk {ky) = Y\ i e tk ^P(k x ,ky)^j: k , ^ W ~ ^ P (*£ , ky ) = 

= E*.,*i T,j ^e^- k ^e^P(k Xl k y )e-< k P(k' x ,k y ) = £ fe ^, ^5 K . k J k ^P{k x ,k v )e-< k P{k' x ,k y ) = 

= y £ fcte e ik ^-^P(k x ,k y )P(k x , k y ) = A £ fcx e ik ^-^P(k x , k y ) = C lk (k y ). (Al) 

I 



Note that we have used the notation P(k x , k y ) to repre- instead of P k to make the momentum dependence easier 
sent the k-dependent projector onto the occupied states 
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to see. 



Appendix B: Entanglement eigenvalues for two 
occupied bands and four sites 



We then see that we have 4 robust modes at exactly 
1/2 in the entanglement spectrum, or exactly twice the 
difference of negative inversion eigenvalues at the two 
inversion symmetric momenta. 



We will analyze only one case of inversion eigenval- 
ues, i.e. when both inversion eigenvalues at k — dif- 
fer from the inversion eigenvalues at k — tt. For sim- 
plicity and without loss of generality, we particularize to 
Ci(0) = C 2 (0) = 1, CiM = C2M = -1, where the orthog- 
onality relations between the wavefunctions at the same 
k and at different inversion symmetric fc's hold due to the 
opposite inversion eigenvalues. Let the wavefunctions of 
the two occupied bands be il>\{k) and ip2(k). Similar to 
the above the arguments in Sec |II C 3} the entanglement 
wavefunctions take the form (i/i^jm?^) where ipA di- 
agonalizes the operator C + mC{P = |[(1 + m)(P + 
P n ) + p* -p(p + im) + VPk(V - im)}. We expand the 
wavefunction ipA- 

ip A = arVi(O) + 02^2(0) + &1V1O) + MbM- (Bl) 

With this expansion, we have to look at the solutions ipA 
which are in the nullspace of PsiV(V + im) + VPk (V — 
im). Once we have found such solutions, we know they 
have 1/2 eigenvalues for m = 1 since (C + C{P)^a = 
\{Pq +P- n )4>A = \^>A due to the fact that wavefunctions 
at different inversion symmetric momenta are orthogonal 
if their inversion eigenvalues are different. We denote the 
overlaps: 



(^i(f)|^i(0)) = a 1 , 

(M?)\M*)) = <*a, 
(V 2 (f)|V'i(0)>=A, 
<V> 2 (f)|V'i(7r)>= / 33, 



(</>l(f )|lM7r)> = OLA 
(iMi)llMO)) = 02 

fe(f)|^W)=A (B2) 



These are the only unknown overlaps, as the wavefunc- 
tion at 37r/2 is related to the one at n/2 by inversion 
symmetry. Thus its overlaps with eigenstates of V are, 
up to a sign, identical to the ones above. We find the 
following two exact half modes of the entanglement spec- 
trum: 

(ai, a 2 , a 3 , a 4 ) = 
(i(/3 3 a 4 - ^403), 0, P4011 - Pioti, /3 3 ai - /3ia 3 ); 
(/3 2 a 3 - /3 3 a 2 , /3 3 ai - Pia*, 1(^1 ~ /?ia 2 ), 0) (B3) 

The other two mid-gap entanglement eigenvalues are ob- 
tained when m = — 1, in which case we have to diagonal- 
ize the operator \[P«V(P - i) + VP%(V + i)}. In this 
case, we expand the eigenstate: 

\M = (ci + c 3 v) \M%)) + ( c 2 + c 4 v) (B4) 

we find 1/2 modes if, just like in the previous section: 

(ci,c 2 ,C3,c 4 ) = (0,1,Q,»), (1,0,1,0) (B5) 



Appendix C: Explicit Proof for the generic two-site 
problem with TV occupied bands 

We now show that the two-site problem with n\ neg- 
ative inversion eigenvalues at k — and n 2 nega- 
tive inversion eigenvalues at k — tt out of a number 
N of occupied bands contains 2|m — n 2 | zero modes 
in the entanglement spectrum. Without loss of gen- 
erality denote the eigenstates of the original Hamilto- 
nian |^i(0)) . . . \ip ni (0)) as the ones with negative in- 
version eigenvalue at k = 0, |-0m+i(O)) . . . \ipjsr(0)) as 
the ones with positive inversion eigenvalue at k = 0, 
1^1(71")) . . . \tpn 2 (j r )) as t ne ones with negative inversion 
eigenvalue at k = tt, |'0ri2+i( 7r )) • ■ • IV'aK 71 ")) as the ones 
with positive inversion eigenvalue at k = tt. Bands at 
the same momentum are all orthogonal, while bands at 
different momenta are orthogonal if they have opposite 
inversion eigenvalues. The simplest case of the above, 
which should be obvious from our previous examples, is 
that of all the N inversion eigenvalues at k = being 
identical and negative of the N eigenvalues at k — tt. In 
this case, the projector at one of the inversion symmetric 
fc's annihilates all the eigenstates at the other inversion 
symmetric k, and the 2N occupied eigenstates of the orig- 
inal two-site Hamiltonian are also the eigenstates of the 
entanglement spectrum at fixed eigenvalue 1/2. Due to 
their orthogonality, they are linearly independent. From 
here it is clear that our formula is physically correct: 
adding the same eigenvalues to both k — 0, tt cannot 
change the result. 

We again expand the eigenstates |-0a) of the correlation 
function Cl = (Pq + Ptt)/2 as a sum over all the occu- 
pied eigenstates, even though these might not be (and in 
general are not) orthogonal: 

V'A = Em=l a m %p m (0) + X)m=ni+1 a m1pm(0) + 

+ Em=l K^mW + Em=n 2 +1 h m^m{^) (Cl) 

In the generic case, we assume that the norms that 
are not fixed to vanish by symmetry (such as differ- 
ent inversion eigenvalues) are all nonzero. In gen- 
eral, it might be the case that not all the eigenstates 
in the expansion above are linearly independent i.e. 
V-Vn (0) might not be linearly independent from a sum 
of the eigenvalues ip m (^) which have identical inver- 
sion eigenvalues. In building the matrix to be diago- 
nalized, we take this into consideration, but when writ- 
ing the eigenvalue equation, we assume they are lin- 
early independent - generically, they will be, because 
there are 2N wavevectors of 2N components. The vec- 
tor (oj, . . . , a ni , a„ 1+ i, . . . , a^v, 6 1; . . . , b„ 2 , 6„ 2 +i, . . . , 6jv) 
has to diagonalize the matrix: 
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2 m x m 
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n,2 x ni 




B 



2 N—nixN—nx 





At 



N-n 2 X JV— ni 



B ni xn 2 



1 

2 n 2 xri 2 







^4 TV— niXjV— n 2 



1 

2 N—niXN—ri; 



where _B, 



1. 



/ 

(C2) 

i = 

.,7V, 



1, . . . , n 2 and A$j = (^(O)|'0i(7r)) with i = + 1, . . 
j ' = ^2 + 1, ■ • • ; N. It is easy to see that this matrix has 
2 1 71,1 — 712 1 eigenvalues at exactly | irrespective of the A, i? 
matrices. We show it for — 0, the generalization to 
n 2 7^ being straightforward. For n 2 = 0, the matrix 
reads: 



/ I 

2 m Xni 



V 





1 

2 N—niXN—n-i 





AiV-TiiXJV 
1 

2JVxJV 



(C3) 



Half of the entanglement eigenvalues at 1/2 are obvious - 
they are the eigenvalues of the ipi(0), i = 1, . . . , n\ eigen- 
states. The remaining eigenvalues must then be part of 
the eigenvalues of the matrix: 



R 



JV— ni 
N 



2N-mxN- ni A N-mxN 
it 1 



.4 



(C4) 



NxN-ri! 



2NxN 



where we have indexed the matrix by the dimension N — 
ni of the upper block-diagonal square matrix and by the 
dimension N of the lower block diagonal square matrix. 
We need to compute the determinant of: 



R 



JV — 711 

N 



(o — A)iv-nixJV- 



.4 



.4 



t 

jvxJV-m 



N-mxN i rnK\ 
! n (C5) 
2 — AJ/VxJV 



We can prove that this matrix has n\ eigenvalues inde- 
pendent of what the matrices A are: as such, we denote 
by Mg a matrix of the form above, but with any ran- 
dom numbers instead of the matrix made out of norms 
matrix Aij — (ipi(0)\i/jj(Tr)). We want to compute the 
determinant of M N ~ ni . By expanding first on the last 
column of the matrix, then immediately after, expanding 
on the last row of all the matrices obtained, we find the 
recurrence relation: 

dct(M^- ni ) = (i - A) clct M^Zx 1 + x • dct(A/^I 1 " 1 - 1 ) 

which, applied successively, leads to: 

det(M^- ni ) = (|- A) r det M^ 1 + 

+ El=i **■(§- A) 1 " 1 det(M^r _1 ) (07) 

where r < N — 1 and numbers. We choose r = 

N — 1 Notice we have so far applied the recursion relation 
only to the lower index of the matrix. We now apply it 



to the det(Af 



JV— 711- 

JV-i 



to obtain: 



det(M^-" 1 - 1 ) = (| - X)P- i detM^~ 1 + 

+ E£ * v ■ (| - A)^ 1 det(M^- 2 ) (C8) 



where p > i is an integer and the x q are complex numbers 
not necessarily equal to the cc,'s in the previous recursion. 
We again choose p = N — 1, which means i < N — 2 in 
the above formula. By plugging in det(M^I™ 1 - 1 ) into 



Eq C7 , and separating the sum to take into account the 



restriction i < N — 2 for which Eq C8 is valid, we find 



det(i<-" 1 ) = 0((i-A) Ar - 2 ) + 



EJV-l ypJV-l 
<70 = 1 "91=1 



JV-1-® 



u qo Uj qi 



(i-s, 



90.JV-1 



: det(Afi 



) 



>2 'V """^"JV-^o+qiV 

where 0((^ — A)^ -2 ) means terms proportional to at 
least a the ( | — A) 2 . By applying the recursion relation 
successively, we obtain: 

det(M^ ni ) = 0{{\ - \) N -^) + 

y^JV-l y^JV-l-go v^JV-l-(g +gi) ^N-l-(q +.-+qi-i) 

(1 - 5 go ,jv-i)(l - <W<?i,JV-i)— (1 - 6g 0+ ... +gi _ uN - 1 ) ■ 
. (i _ A) (*+..-h B -(i+i) defcCM^;?^ j) (C9) 



when iV — ni — (Z + 1) = 0, the matrix M 

fully diagonal and has determinant (i — A) 
Hence the term in the sum is: 



JV-m-(z+i) 

N-(q + ...+qi) 



IS 



JV-( 9o +.-+gi) 



(h—X) qo+ — +qi ~ l ~ 1 (~—\) N ~' IO ~ 



"« = (--Ar (CIO) 
\JV-(J+l)^| 



while the lowest order term in 0((i — A) 7V is 
again (| - A)™ 1 . We hence see that the determinant of 
det(M^~ ni ) ~ (| — A)" 1 , and these are the remaining m 
eigenvalues. Combined with the first n\ eigenvalues, we 
see we have a total of 2n\ eigenvalues in the entanglement 
spectrum at 1/2, just as the formula predicts. 



Appendix D: Properties of the ID adiabatic 
connection under inversion symmetry 

We begin by recalling some basic facts of inversion 
symmetric Hamiltonians. Assume \m{k)) is an eigenstate 
of the Hamiltonian at energy 

H(k)\u i (k)) = E i (k)\u i (k)) (Dl) 

Then V \ui(k)) is necessarily an eigenstate at — k of the 
same energy: 

VH{k)V~ l V \u t (k)) = H(-k)V \ Ul (k)) = VH(k) \ui(k)) 
= PE i (k)\u i {k))=E i (k)P\u i (k)}. 

We assumed no degeneracies in the spectrum which 
means 



M-fc)) =e^V\ Ul {k)). 



Thus we have: 



\ Ul (k)}=e- ia *V\u m (-k)) 



(D2) 



(D3) 



with Ei(k) = Ei(-k). 

We assume we are in an insulating state where the 
charge polarization is 



3G 




Using Eq. |D3| we get 

p l = i E Je occ. {Vui{-k)\ e^d k e-^V \ Ui (-k)) = 

= -iej^ f E. ieo cc. CPuii-k)] d k V M-fc)) - ief^ f E ieo cc.H)^- (D5) 

The last term is an integer because it depends only on a, — a_ 7r which can at most be 2nj, so we will drop it and be 
left with: 

Pi = i E ieo cc. (ui(-k)\d k \ Ui {-k)) = ie /"* * £ ieo cc. ( Ui (fc)| d_ fc | Uj (fc)> = 

= -*e/! w f E ieo cc. 9- k \ui(k)) = ief^ f £ l6 occ. 9 fc | Ui (fc)) = -P x . (D6) 

Since Pi is defined only mod m the two values consistent with inversion symmetry are: 

A=0,|. (D7) 

We define the trivial insulator as Pi = 0, while the topological one is Pi = |. The e/2 value means that when the 
system is cut in two, there exists half a charge on each end. 



Appendix E: Relation between adiabatic connection 
and B(k) in ID 

We now prove that the matrix 

Bij = {ui- k \V\uj tk ). 

is unitary and that 

i(-fc) = -A{k) + iTr[B(k)V k B^(k)] 

Wc define the matrix B^ as the matrix connecting the 
bands at k with the ones at —A;: 

\u i (-k))=B* j (k)V\u j (k)) (Ef) 

where i,j run over the occupied bands 1, ...,N. From 
this formula we see that 

That Bij has to be unitary can be easily seen by ex- 
tending the matrix to belong to all bands, occupied and 



I 

unoccupied. Since we have a full gap in the system none 
of the occupied bands at k can transform to unoccu- 
pied bands at — k and vice- versa (otherwise we would 
not have an insulator). This means that the full ma- 
trix P fu11 = diag(B occu P' lcd , B unoccu P icd ) is block diagonal. 
One can prove the full B is unitary by using complete- 
ness: 

J2 \u. l (k)}(u l (k)\=I (E2) 

i£all bands 

we have (double index means summation) 

{B^B)ij = (B*) im B mj = (u i:k \V^ \u m ,-k) (u m ,-k\'P\uj,k) 
= (u^klV^V \u j:k ) = (u i:k \ \u jtk ) = Sij 

Hence, since the full matrix B is unitary, so are the 

(un)occupicd 

We now want to express a connection between A(— k) 
and A(k). We have 



i(-fc) = -i {u it - k \V- k \ui,- k ) = i{ui(-k))* a V k (ui(-k)) a = iBuityV^iui^VkiB^ifyVaeiUj^e) = 

= iB u (k)\ykB^(k)](ui )k )*p5p0(Uj,k)e + iBii(k)B^{^ = 
= iBuWlVkBfjikWn + iBu(k)B^(k) {u hk )}V k {u jtk ) = 
= iB u (k)\^ k B* a {k)] +id 3 i {ui, k )* p V k (uj. k ) p = iBa(k)\VkB^(k)] +i (uj, k )* V k (uj !k )p = 

= iTr[B(k)V k B\k)] - A(k) (E3) 

where repeated indices are summed over. 
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Appendix F: Magneto Electric Polarization as where the non-Abclian adiabatic connection is is 

Winding Number of the B(k) matrix A^ (k) = — i(u a>k \V ki \up jk ). The above implies that 

The Abelian A(k) in the previous section obeys a spe- 
cial case of the more general non- Abelian transformation 



A(-fc) = -PA(fc)P t + iB{k)VB\k) (Ff) 
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the Berry gauge fields at k and —k are non-abelian gauge 
transformed of each other. The field strength gauge 
transformation is: 



F ij (-k) = B(k)F ij {k)B^(k) 



(F2) 



From here the magneto-electric polarizability is easy but 
tedious to obtain: 



Pa 



l 

16tt 2 



/ d 3 fce iifc Tr[(4(fc) - liMk)A 3 (k))A k (k)] 



1 

16tT 2 



/ (fkdjkTrKFiji-k) 



i 

16tt 



1 

167T 2 



Jd 3 k 



jiM-k)A 3 (-k))A k (-k)] = 
J d 3 ke ijk Tr[(B(k)Fij(k)B^ - |i(S(fc)i i (fc)i J -(fc) J Bt(fc) - iB{k)A l {k)d 3 B^ (k) - 
+ i{d i B{k))A j {k)B^ (k) + d l Bd J B^){-B{k)A k {k)B^{k) + iB{k)d k B^ \k))\ = 
e l]k Tr[(F l3 (k) - liAii^Ajik))!^)} - ^ f < Pke ijk Tr[(B(k)d i Bl)(B(k)d j Bl)(B(k)d k Bt)] 



which proves the formula in the text. 



Appendix G: Proof of Fu-Kane Formula 3 

In this appendix we provide an alternative proof for 
the Fu-Kane formula for the Z 2 invariant of 3D T and 
V invariant insulators. Consider the Bloch Hamiltonian 
H (k) of an insulator with both inversion V and time- 
reversal symmetry T. We have 



1. 



T 2 



TH{k)T- 1 
VH{k)V- 1 



-1, [P,T] 

H(-k) 

H(-k) 







and hence: 



VTHityiTT)- 1 = H(k), (VT) 2 = -1 



(Gl) 



(G2) 



which proves that the insulator has doubly degenerate 
bands at each momentum k. This does not depend on 
the dimensionality of the space. The two ingredients we 
will use to prove Eq. 108 are (i) band crossing arguments 
between a 3D insulator in the trivial atomic limit and a 
topologically non-trivial insulator, and (ii) the fact that 
a 4 + 1-d Dirac Hamiltonian changes its 4D 'Hall con- 
ductance' by 1 if there is a band crossing of four bands 
- these are actually two doubly degenerate bands. With- 
out loss of generality we consider a topological insulator 
with four bands - two occupied bands and two unoccu- 
pied. As we take a generic T and V invariant insulator 
through a phase transition, four ( i.e. two doubly degen- 
erate) bands are generically needed. 

Assume we have two T and V symmetric Hamiltonians 
in 3D hi(k) and h 2 (k) with: 

vkM^r- 1 = kA-k). 

We choose h\(k) to be trivial (in the atomic limit with 
all hoppings taken to vanish). We now construct a time- 
reversal and inversion invariant interpolation between 
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these two Hamiltonians. We first prove that such a 
gapped interpolation exists. To see this, it is easiest to re- 
main with our reduced Hamiltonian with 4-bands, which 
represents the generic effective Hamiltonian of the two 
doubly-degenerate bands immediately above and below 
the Fermi level, out of the total N bands in the insu- 
lator. In this basis, the 3D effective insulating Hamil- 
tonian with inversion and time-reversal symmetry has 
co-dimension 2 (there are 3 momenta k x ,k v ,k z and 5 
Clifford generator T a matrices in which a 4-band insu- 
lating Hamiltonian with doubly degenerate bands can be 
expanded). In AD, a topological insulator with inversion 
and time-reversal symmetry has co-dimension 1 (there are 
4 momenta and 5 T a matrices), so there is still always a 
way to make it gapped. This shows that generically, a 
3D or AD insulator with inversion and time-reversal is al- 
ways gapped. It can be made gapless by tuning 2 and 1 
parameter(s) respectively. Now, let the gapped interpo- 
lation between h\(k) and h 2 (k) be h(k, 9) which satisfies 
the properties: 



h{k,0) = h^k) ; h(k,ir) 
Thfce)^ 1 = h(-k,-6 
Vh{k,ff)V~ x = h(-k,-6 



h 2 (k) 



(G3) 



The interpolation between the 3D Hamiltonians is cho- 
sen this way so that if one interprets 9 as a fourth mo- 
mentum then the resulting AD Hamiltonian would re- 
spect inversion and time-reversal. It was shown in Ref. ,9 
that 2(P 3 (h 2 {kj) - P 3 (hi(k))) = C 2 (h(k,9)) where C 2 is 
the second Chern number of the AD Hamiltonian h(k, 9). 
Since we chose h\{k) to be trivial 



P 3 (/ii) = mod n e Z. 



(G4) 



Hence, if the second Chern number of the 4-D h(k,9) is 
odd, we have 



P 3 (ft 2 ) = l/2 modneZ 



(G5) 
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giving rise to the result that if C 2 of the AD Hamiltonian 
is odd, then either hi(k) or h 2 (k) is a non-trivial topo- 
logical insulator. Since we pick hi(k) to be our reference 
trivial Hamiltonian then this would imply that /12(A) is 
non-trivial. 

To understand how to classify the 3D insulators we 
first need to understand how to get a AD insulator with 
an odd second Chern number. The 4D trivial Hamil- 
tonian is simply a momentum-independent interpolation 
between h\(k) and itself. This clearly has vanishing C 2 . 
Since C2 is a topological invariant, we must have a gap- 
closing phase transition to change it. As the system is 
inversion and TR invariant, we have to analyze the cross- 
ings between two doubly-degenerate bands. This is the 
generic case, even for insulators with an arbitrary number 
of bands N, because we can build an "effective" Hamilto- 
nian close to the transition which will be a 4-band model. 
With time-reversal and inversion, the Bloch Hamiltonian 
has to be of the form H(k) = d a (k)r a . We first con- 
sider transitions which occur away from the invariant 
momenta. Because of inversion (or time-reversal) a gap 
closing at k must be matched by one at —k. Such tran- 
sitions can be tuned by a single parameter. Since the 
gap closing and re-opening happens away from an in- 
variant momentum, the inversion eigenvalues of the oc- 
cupied bands remain unmodified. The non-abelian adia- 
batic field strengths at the two k-points are equal (up to 
a gauge transformation, and a minus sign in case of time- 
reversal). Thus, a gap closing at two points, k and — k 
makes the total change in second Chern number which is 
even. Thus a 4D Chern insulator with inversion and time 
reversal has even second Chern number if the inversion 
eigenvalues of the occupied bands are the same as in the 
atomic limit. 

We now look at the case where the gap closing and re- 
opening happens at an inversion symmetric point. In this 
case, the Hamiltonian is still 4x4 but the co-dimension 
of any two doubly degenerate bands is 5 because we have 
5 Clifford matrices and no tunable momenta - the mo- 
menta are fixed at the inversion symmetric points. If 
the bands involved in the phase transition at the inver- 
sion symmetric points have the same inversion eigenval- 
ues (e.g. all positive) then, up to a gauge choice, the 
inversion matrix is the identity operator. This means 
that [P,H(0)] is trivially satisfied and does not provide 
an additional constraint and so the crossing is always 
avoided with such a large co-dimension. However, if the 
bands have different inversion eigenvalues, (e. g. bonding 
and anti-bonding bands), so that the inversion matrix in 
a specific choice of basis is V = 12x2 ® i~ z , then we find 
that the effective Hamiltonian matrix at the inversion 
symmetric points Gj/2 is 



H(Gi/2) = M{P 



(G6) 



second Chern number will change by ±1 if only one mass 
goes through zero, or more generally, it will change by 
the number of masses that go through zero (multiplied 
by a sign) and (ii) the inversion eigenvalues of a pair 
of occupied bands change sign as a result of every Mj 
that switches sign. This is the crucial difference be- 
tween phase transitions which change the first and sec- 
ond Chern numbers. In two dimensions a change in the 
first Chern number is accompanied by a single inversion 
eigenvalue switch since the minimal crossing is between 
two-bands. In 4D, changes in the second Chern num- 
ber require a four-band crossing and thus two inversion 
eigenvalues are exchanged. This carries over to 3D where 
we are interested in the Z 2 valued parity of the second 
Chern number. A change of the parity is effected when a 
4-band crossing switches two inversion eigenvalues. We 
hence proved that a 4D Chern insulator with inversion 
and time-reversal symmetry has odd second Chern num- 
ber if the product of half of the inversion eigenvalues (half 
meaning each Kramers' pair is only counted once) at all 
TR invariant points is —1. By keeping h\(k) constant and 
changing only h 2 (k) through a gap closing and reopening 
we can make h(k, 9) a 4D insulator with an odd second 
Chern number. Thus h 2 (k) is a non-trivial Z 2 insulator 



when Eq. 108 is negative and we have proved the Fu and 
Kane formula. 



Appendix H: Proof that the winding number is 
independent of a smooth U(l) phase 

We begin with the generic form of a U (2) block of the 
entire B(k) matrix 

B(k) = e i ^(f(k)I + ig a (k)a a ) (HI) 



where 



(f(k)) 2 +g a (k)g a (k) = l. 



(H2) 



where Mi is the mass at the inversion symmetric point 
ElHH] As such, a gap closing transition whereby the 
Mi changes sign is accompanied by two effects: (i) the 



We drop the argument (k) which is implied from now on. 
By using: 

O a Ob = Sab + KabcVc (H3) 

we find: 

B^diB = i [d^ + (fdig c - g c dif + e bac g b dig a )a c ] (H4) 
Define 

G ic = {fd t g c - g c dif + ebac9bdig a ) (H5) 

so that: 

B%B = 1 [di4> + G ic a c ] . (H6) 
The winding number density is: 
n w = e ljk Tr[B^ diBB^ djBB^ d k B] 

= -i^ijkTr [(di<p + G ia a a )(dj(f) + G jb a b )(d k (f> + G kc a c 

Thus the U(l) phase does not contribute if it is smooth. 
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Appendix I: Proof that the 17(1) phase is smooth in 

3D 

We have seen that the B(k) matrix decomposes into 
trivial one dimensional and non-trivial 2 dimensional di- 
agonal blocks. Each 2 dimensional block can be seen as 
map from T 3 into U(2). Here we show that the U(l) 
phase of these maps is smooth and for this reason it can 
be taken out from winding number calculation. For this, 
we will show that a discontinuity in the phase implies 
the appearance of the (_) eigenvalue pattern of inver- 
sion eigenvalues at one fc inv point. Since this pattern is 
excluded for the U(2) blocks that cannot be reduced to 
U(l) x £7(1), it follows that the phase is actually smooth. 

Now, a simple example of a map with a non-smooth 
phase, to keep in mind for intuition is: 

B{k) = e2 lfc (cos| -Msinf a z ) = e ie ^U(k) (II) 

where 17 (ft) G SU{2). Note that this B{k) satisfies all 
the required constraints and we can see that in Eq. [IT] 
everything is smooth as k is continuously varied, except 
when we are close to the end of the 2ir cycle, that is when 
we close the loop around the Brillouin torus. If we start 
for example from k=— ir, one can see that the phase and 
the SU(2) factor take different values when k approaches 
7r, but as a whole B{k) reaches the same value as the one 
we started with at fc=7r. It is important to note that 
the discontinuity in the phase cannot be removed by any 
continuous deformation of B(k). 

In the following we argue that such a discontinuity in 
the phase immediately implies the existence of a (_) pat- 
tern of inversion eigenvalues at some ft inv point. As we 
already pointed out, the non-removable discontinuities in 
the phase are entirely due to the nontrivial topology of 
the Brillouin torus. The question is what happens when 
k is taken over a nontrivial loop of the torus? So let us 
start from a k} av point and continuously vary k until it 
advances by 2tt along one of the ki, k 2 or fc 3 directions, 
that is, we are taking k around non-trivial loop of the 
torus. As before, we will work with the following expres- 
sion of the U (2) matrix: 

B{k) = e i8h (fk + igko). (12) 

At the fc inv points, 9 can take only integer or half-integer 
values (in units of 7r), and an important observation is 
that, if 6 is half-integer at a fc inv point, then B(k inv ) nec- 
essarily takes the form ha (n= unit vector), which has 
±1 inversion eigenvalues. So let us assume that 8=0 at 
kf nv (otherwise the (1) pattern will show up). After we 
complete the loop around the torus, in general, we can 
have: 

0*£, r +2* = + { (2nTl)n (b) (I3) 

In case (a) the phase factor returns to its original value so 
it is smooth, in which case there is no problem removing 



it from the winding number calculation. In case (b), 
the phase factor changes sign. Now let us consider the 
inversion properties of B(k) relative to the kf av =k} nv +Tr 
point, which also belongs to the loop. We have: 

e 6k ?n V + k (f k 2 +k + ig k 2 +k a) 
-id ^ 

Since f k and g k take only real values, we can see that if 
6 k 2 +k is half-integer then necessarily 9 k 2 _ k is also half- 
integer. In other words, the half-integer values of 9 k come 
in pairs, unless k is the fcr^ invariant point itself. Since 9 
varies from to (2n+l)7r as we vary k by 27r, there will be 
an odd number of times when 9 k assumes a half-integer 
value. It is then evident that one such half-integer value 
must occur at kf nv and consequently the ( — h) inversion 
eigenvalue pattern will show up at kf m . 

Appendix J: Inversion eigenvalue patterns and the 
corresponding magneto-electric polarizability 

In this appendix we provide some explicit examples of 
our proof of the connection between inversion eigenval- 
ues and the magneto-electric polarization for inversion 
invariant insulators. We focus on 4-band, gapped Hamil- 
tonians H(k) in 3-dimensions displaying various patterns 
of inversion eigenvalues. Assuming two occupied bands, 
we will write these inversion eigenvalue patterns as 

+ + + + + + ++ ( m 
+ + ++, lJij 

for example, where one should understand that the two 
non-zero eigenvalues of Pk inv V Pk iIlv are +1 and —1 at 
four fc inv points and +1 and +1 at the remaining ft inv 
points. The order of the k inv is not important for the 
magneto-electric polarizability, but to be precise we will 
order them as 

{ (0, 0, 0), (tt, 0, 0), (0, tt, 0), (0, 0, tt), (tt, tt, 0), (tt, 0, tt), 

(0,7T,7t), (7r,7T,7r)}. 

For each case, o compute the magneto-electric polariz- 
ability P 3 , we will begin with a trivial reference Hamil- 
tonian (P3 = 0) and find a gapped and inversion sym- 
metric interpolation H(k, 9) between this and our ex- 
ample Hamiltonians. We can compute P3 as half the 
second Chern number C2 generated by the interpolation 
H (ft, 9) which is considered as a 4+1-d inversion symmet- 
ric Hamiltonian. An odd C2 corresponds to a non-trivial 
insulator with half-integer magneto-electric polarizabil- 
ity, while an even C2 indicates a trivial insulator. This 
method of calculation is convenient because Ci can be 
computed using a gauge-invariant projector method , so 
we can bypass the task of finding a smooth gauge. All 
of our numerical observations are in agreement with our 
mathematical proofs above. 
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1.) The pattern 



-+ 
-+ 



(J2) 



is seen with the following gapped Hamiltonian: 

H{k) = sin fciTi + sin k 2 T 2 + sin k 3 T 3 

+(— 2.5 + cos ki + cos fc 2 + cos k 3 )To (J3) 

+0.75(cos fci + cosfe 2 )(f 15 + f 25 - f 35 - f 12). 

This Hamiltonian can be connected through a gapped 
interpolation with the Hamiltonian: 

Ho(k) — sin kiTi + sin k 2 Y 2 

+ (-2.5 + cosfci + cosfc 2 )f (J4) 

+0.75(COS fci + COSfe 2 )(r 15 + f 25 - f 35 - f 12 ). 

using the inversion symmetric homotopy 
H(k,9) = |(1 + cos9)H(k) + §(1 - cos9)H {k). (J5) 

The second Chcrn number generated by H(k, 9) is C 2 =0 
(even) and the P 3 of H (k) is zero because it depends only 
on k\ and k 2 . Consequently, the P 3 of H (k) is mod(Z). 
2.) The pattern 



-+ 
-+ 



(J6) 



is seen using the following gapped Hamiltonian: 

-ff (fc) = sin feiTi + sin fc^I^ + sin k 3 T 3 
+(—2.5 + cos ki + cos k 2 + cos k 3 )T 
-(2 - cos fc 3 ) (cos ki + cosfc 2 )(f 15 + f 25 - f 35 - fi 2 ). 

This Hamiltonian can be connected through a gapped 
interpolation with the Hamiltonian: 

Ho(k) = sin fciTi + sin fc 2 T2 
+ (-2.5 + cos fci + cosA;2)fo (J8) 
-2(cosfci +cosfc 2 )(fi 5 +f 25 -f 35 -fia). 
using the inversion symmetric homotopy 
H(k,9) = \{l + cos9)H(k) + |(1 - cos9)H {k). (J9) 

The second Chcrn number generated by H(k, 9) is C 2 =Q 
(even) and the P 3 of Ho(k) is zero because it depends only 
on fci and k 2 . Consequently the P 3 of H{k) is mod(Z). 
3.) The pattern 

----1111 («o) 



is seen with the following gapped Hamiltonian: 
H (fc) = sin fciTi + sin fc2T2 + sin fc 3 T3 



+(— \ + cos fci + cosfc 2 + cosfc 3 )f + |f 



(Jll) 



25- 



This Hamiltonian can be connected through a gapped 
interpolation with the trivial Hamiltonian: 

Ho{k) = smfcifi +sinfc 2 r 2 +si n fc 3 r3 — k^o, (J12) 



using the inversion symmetric homotopy of Eq. |144| 
The second Chern number generated by -ff (fc, 9) is C2=2 
(even) and consequently P 3 = mod(Z). 
4.) The pattern 



+ + + ++ 



(J13) 



is seen using the following gapped Hamiltonian: 

H (fc) = sin fciTi + sin fc2l2 + sin fc3T3 — |fo 

+ (2 + cos fci + cosfc 2 + cosfc 3 )f 15 (J14) 

+ ( — 2 + COS fci + COS fc 2 + COS fc 3 )f25 

This Hamiltonian can be connected through a gapped 
interpolation with the trivial Hamiltonian: 

H (k) = sinfcifi + smfc 2 f 2 + sinfc 3 f 3 - |f (J15) 

using the inversion symmetric homotopy 

H (fc, 9) = sin fciTi + sin /^I^ + sin k 3 T 3 — |lo 



+ |(1 + COS0)(2 + COS fci + COSfc2 + COS fc3)f 15 
+ 5(1 + C0S9)( — 2 + COsfci + COsfc2 + COS fc3)T25 

+ sin^ f 5 . 



(J16) 



The second Chern number generated by H(k, 9) is C 2 —0 
(even) and consequently P 3 = mod(Z). 
5.) The pattern 



h + + + ++ 

h + + + ++ 

is seen using the following gapped Hamiltonian: 



(J17) 



H (fc) = sin fcif 1 + sin k 2 t 2 
+(-§ + cos fci + cosfc 2 )f (J18) 
+(-| +0.1(cosfci -cosfc 2 ))f 25 . 

This Hamiltonian can be connected through a gapped 
interpolation with the trivial Hamiltonian: 

H (k) = sinfcifi + sin fc 2 f 2 + sin fc 3 f 3 - |f (J19) 
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using the inversion symmetric homotopy 

H{k) = smfciTi + sinAfef^ + |(1 — cosO) sinA^r^ 
+(-§ + |(1 + cos6»)(cosfci + cosfc 2 ))f 
+§(l + cos0)(-| +0.1(cosfci -cosfc 2 ))f 25 
+ sin# r 5 . 

(J20) 

The second Chcrn number generated by H(k, 8) is 6*2=0 
(even) and consequently P3 = mod(Z). 



Appendix K: Response Theory Argument in Higher 
Dimensions 



Although we do not rigorously prove anything in higher 
dimensions we note that the physical response arguments 
presented in this work and Refs. [3] and [55] continue to ap- 
ply in higher dimensions. This gives us a hint that there 
are interesting inversion symmetric topological insulators 
in higher dimensions. For even spacetime dimensions the 
topological response actions all take a standard form 



S eff[ A u 



d 2n xP n e° 



.F„ 



where P n is a response coefficient and is the elec- 
tromagnetic field-strength tensor. If these were all dy- 
namical fields then the entire action must transform like 
a scalar which means that the intrinsic response coeffi- 
cients P n must transform the same way as the product of 
the electromagnetic fields. Thus, P n must be odd under 



inversion symmetry for all n. Additionally for even (odd) 
n, P n is odd under T (C) symmetry. Since we are talking 
about external electromagnetic fields we do not transform 
them under the symmetry operation and hence only the 
P n are changed. Thus in alternating even spacetime di- 
mensions topological insulators are protected by either 
T or C symmetry. However, in every even spacetime 
dimension there is a topological insulator protected by 
inversion symmetry which has a topological response. In 
all of these cases P n is not gauge invariant under transfor- 
mations of the occupied wavefunctions and the oddness 
under the different symmetries quantizes P n to take only 
two independent values. This yields a Z2 classification 
for the response coefficient in line with the arguments 
presented in Ref. HI 

In odd space-time dimensions the generalized Chern- 
Simons terms describing the electro-magnetic responses 
in topological insulators are compatible with inversion 
symmetry. The general action is 



niodd) 
^eff 



(n + l)!(27r) r 



-1 flifl2 - 



■ a 2 n+i A f 1 



■ F 



The quantity C n is the n-th Chern number and is even 
under inversion symmetry and thus not restricted by the 
requirement of preserving inversion symmetry. If n is 
even (odd) then the action is also compatible with T 
(C) symmetry. We have not proven it, but our intuition 
suggests that the parity of C n can be characterized by the 
inversion topological invariants given by the x-p defined 
above when calculated in (2n + 1) — d. 
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